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SUCTION  I 


INTRODUCTION 

Reflections  from  ionospheric  irregularities  can  he  an  important  source  of  clutter  for 
an  IIF  backscatter  radar  located  in  the  midlatitude  with  the  antenna  beam  oriented  towards 
the  polar  region.  Although  theoretical  and  experimental  investigations  of  radar  reflections 
from  K-laycr  auroral  and  I'-laycr  ionization  irregularities  have  been  conducted  extensively 
at  frequencies  between  the  IIF  and  U1IF  band,  only  limited  experimental  data  on  the  statis- 
tical characteristics  of  IIF  ionospheric  echoes  are  available  which  would  be  directly  appli- 
cable for  IIF  backscatter  radar  design  and  performance  determination.  However,  sufficient 
information  and  knowledge  concerning  ionospheric  backscatter  reflections  do  exist  to  pro- 
vide, to  a first  approximation,  a model  for  the  synthesis  of  an  HF  backscatter  radar. 

Radar  reflections  from  auroral  ionization  occur  in  a region  which  satisfies  the  condi- 
tion of  near  orthogonality  of  the  incident  radiation  with  the  magnetic  field  lines  of  force, 
the  amplitude  of  the  auroral  echo  being  aspect -angle  sensitive  (Leadabrand,  et  al. , 1967; 
Chesnut,  et  al.  , 1968;  Hates  and  Albcc,  1969). 

Although  radar  returns  from  auroral  ionization  at  VHF  and  UHF  arc  normally  confined 
to  E-layer  heights  of  about  110  km  (Unwin,  1959;  Watkins  and  Sutcliffe,  1965;  Leadabrand, 
et  al. , 1965;  Abel  and  Newell,  1969),  106-MIIz  auroral  echoes  have  been  detected  as  high  as 
300  km  (Schlobohm,  et  al.  , 1959). 

In  the  case  of  IIF  radar  transmissions,  reflections  extending  up  to  F-layer  heights, 
i.  e. , 300-400  km  altitude,  have  been  observed  as  well  as  in  the  E-layer  (Malik  and  Aarons, 
1964;  Bates  and  Albec,  1969,  1970;  Baggaley,  1970;  Au  and  Hower,  1970). 

In  addition  to  the  F-layer  reflections,  an  I1F  radar  should  be  capable  of  detecting 
ionospheric  echoes  from  over  the  horizon  as  a result  of  the  normal  ionospheric  refraction 
phenomenon.  The  HF  radar  backscatter  soundings  by  Bates  (1969,  1970)  in  Alaska  con- 
firmed the  presence  of  over-the-horizon  ionospheric  echoes  at  extremely  long  ranges  cor- 
responding to  regions  in  space  where  the  propagation  direction-magnetic  field  line  perpen- 
dicularity condition  is  satisfied. 
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A number  of  theories,  which  are  summarized  by  Hooker  (1900)  and  Moorcroft  (19'Jla, 
19011))  have  evolved  invoking  magnetic  field-aligned  Ionization  to  explain  the  aspect  sensi- 
tivity. The  earliest  theories  of  the  HF  radar  aurora  suggest  field  aligned  sheets  of  ioniza- 
tion with  critical  frequency  exceeding  that  of  the  transmitted  frequency.  Subsequent  obser- 
vations at  VHF  and  UHF  resulted  in  apparent  critical  frequencies  much  higher  than  physically 
realizable.  Booker,  et  al.  (1955)  postulated  that  the  auroral  reflections  occurred  from 
columns  of  ionization  with  critical  frequency  below  the  operating  frequency  (weak  scattering). 
The  investigation  by  Lyon  and  Forsyth  (1902)  indicated  that  the  critical  frequencies  in  the 
auroral  ionization  columns  may  be  as  high  as  100  MHz. 

The  analysis  of  radar-auroral  clutter  data  collected  at  six  frequencies  from  50  to 
3000  MHz  simultaneously  at  Homer,  Alaska,  by  (’hesnut,  et  al.  (1908)  has  revealed  the 
following: 

1.  Aspect  sensitivity  is  nearly  independent  of  frequency; 

2.  Scattering  cross  section  decreases  about  10  dB  per  degree  of  aspect  angle; 

3.  Scattering  cross  section  varies  exponentially  with  frequency,  the  slope  being 
time-variant  with  an  average  of  33  dB  per  1000  MHz;  and 

4.  Thcro  is  a moderate  correlation  between  the  radar  aurora  and  the  visual  aurora. 

The  existence  of  two  types  of  auroral  echoes,  which  have  been  termed  diffuse  and 
discrete,  have  been  reported  (Prcsnell  et  al.  , 1959;  Leadabrand,  1960).  The  discrete  echoes 
appear  to  result  from  auroral  ionization  which  is  restricted  in  space  to  a small  reflecting 
region  and  is  prevalent  at  night.  Diffuse  echoes,  on  the  other  hand,  appear  to  originate  from 
large  reflecting  regions  and,  most  often,  are  present  during  the  day.  Discrete  echoes 
exhibit  little  pulse  lengthening  while,  for  the  diffuse  eehoes,  the  pulse  lengthening  could  be 
on  the  order  of  hundreds  of  kilometers. 

In  this  report,  an  estimate  is  made  of  the  characteristics  of  HF  radar  ionospheric 
echoes  reflected  from  the  polar  ionosphere.  The  experimental  radar  data  acquired  at  fre- 
quencies between  HF  and  UHF  are  employed  in  the  analysis.  The  topics  whieh  are  discussed 
are  the  spatial-geographical  extent  of  the  HF  echoes  for  a midlatitude  location,  the  amplitude 
and  backscatter  cross-sectional  area,  the  Doppler  frequency  variation,  the  frequency  of  oc- 
currence, the  diurnal  and  seasonal  variation  and  the  correlation  with  solar -geophysical 
conditions. 
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SECTION  II 

IONOSPHERIC  BACKSCATTER  ECHO  REGION'S 

2.  1 INTRODUCTION 

E -layer  auroral  and  F-layer  echoes  are  the  result  oLradlo  wave  scattering  from 
electron  density  irregularities  aligned  along  the  magnetic  field  lines.  The  characteristics 
of  the  field-aligned  scattercrs  are  such  that  the  radar  echoes  are  observed  from  a small 
range  of  angles  about  perpendicular  incidence  to  the  geomagnetic  field  lines  and  that  the 
echo  magnitude  is  aspect  angle  sensitive.  In  addition  to  the  perpendicularity  requirement, 
it  is  necessary  that  this  geometric  configuration  take  place  at  ionospheric  altitudes,  i.  c. , 

SO  km  and  above. 

The  probability  of  observing  radar  reflections  from  ionization  irregularities  is  also 
dependent  upon  the  frequency  of  occurrence  of  E-layer  auroral,  or  spread-F  irregularity 
activity  in  the  region  of  interest.  Reflections  can  be  expected  to  occur  in  regions  where  both 
the  conditions  of  near-perpendicularity  at  ionospheric  heights  and  high  auroral  or  spread-F 
irregularity  activity  are  satisifed.  The  fulfillment  of  only  one  requirement  is  not  sufficient 
to  warrant  a radar  reflection. 

An  appropriate  model  which  is  used  to  define  the  condition  of  auroral  activity  is  the 
Feldstein  and  Starkov  (1967)  auroral  belt  (oval)  model.  According  to  Bates,  et  al.  (1966, 
1969),  the  spatial  extent  of  the  HF  radio  aurora  compares  closely  with  the  optical  auroral 
ovals  of  Feldstein  and  Starkov. 

Since  the  location  and  extent  of  the  auroral  ovals  are  a function  of  time  and  geomagnetic 
activity  (Feldstein  and  Starkov,  1967)  it  can  be  expected  that  the  auroral  echoes  could  appear 
over  a wide  area  during  magnetically  disturbed  periods  for  an  HF  backscatter  radar  located 
at  high  latitudes. 

The  computational  procedure  discussed  in  Appendix  A is  employed  in  the  determination 
of  the  magnetic  field-aspect  angles  presented  in  this  report. 

2.  2 ELECTRON  DENSITY  MODELS 

The  refractive  index  in  the  ionosphere  is  expressed  by  the  relationship 
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3 -10 

where  N is  the  electron  density  (electrons /cm'  ),  e is  the  electron  charge  (4.  8 x 10  esu), 
e -‘>8 

me  is  the  electron  mass  (9. 1 x 10  gm)  and  w is  the  angular  frequency  of  the  transmitted 
wave  (rnd/s). 

The  distribution  of  electron  density  with  height  is  assumed  to  follow  the  Chapman 
model  of  the  form 
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where  Hg  is  the  scale  height  of  the  neutral  particles  and  Nm  is  the  electron  density  at  the 
level  of  maximum  ionization,  hm.  The  values  of  the  parameters  defining  the  daytime  and 
nighttime  electron  density  profiles  and  the  equivalent  plasma  frequencies  of  the  maximum 
ionization  levels  employed  in  this  analysis  are  presented  in  Table  2-1. 

Plots  of  the  vertical  distribution  of  electron  density  are  shown  in  Figures  2-1  and  2-2. 

For  the  daytime  ionosphere  models,  the  F-layer  plasma  frequency  is  9.3  MHz  and,  for  the 

nighttime,  4.  9 MHz.  The  difference  between  the  two  daytime  profiles  is  the  height  of  the 

peak  of  the  F-layer  maximum.  For  model  A,  h = 250  km  and  model  B,  h = 300  km. 

max  max 

Similarly,  for  the  nighttime  profiles,  for  model  A,  h 250  km  and  model  B,  h =350 

max  max 

km.  For  simplicity  of  data  presentatioi.,  t ic  E-layer  echoes  are  assumed  to  originate  at 
altitudes  below  150  km  altitude  and  the  F-layer  echoes  above  150  km. 

It  should  be  noted  that  the  ionospheric  index  of  refraction  is  also  a function  of  both 
the  electron  collision  frequency  and  the  earth’s  magnetic  field.  In  order  to  simplify  the 
analysis,  both  parameters  are  neglected  in  the  computations. 

For  frequencies  on  the  order  of  10  MHz  and  above,  and  at  altitudes  greater  than  80 km, 
the  effect  of  the  collision  frequency  term  on  the  index  of  refraction  is  negligible  (Davies,  1965). 

The  refractive  index  is  slightly  in  error  when  the  magnetic  field  is  neglected.  The 
maximum  error  in  the  refractive  index  occurs  at  the  peak  of  the  F-layer  and,  at  30  MHz, 
is  less  than  0.4  percent.  At  20  MHz,  the  maximum  deviation  increases  to  slightly  greater 
than  1 percent  (Millman,  1975). 

2.3  E-LAYER  BACKSCATTER  ECHO  REGIONS 

The  contours  of  zero  degree  aspect  angle  (i.  e.  , zero  degree  off-perpendicular  angle, 
propagation  direction  is  orthogonal  to  the  magnetic  field  lines),  fo„  the  E-layer  echoes,  as 
viewed  from  a midlatitude  location,  are  plotted  in  Figures  2-3  through  2-8  for  both  daytime 
ionospheric  models  and  as  a function  of  transmission  frequency  (10,  20  and  30  MHz)  and  the 
various  states  of  magnetic  activity. 
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Figure  2-1.  Chapman  Model  of  Daytime  Electron  Density  Profiles 


Figure  2-2.  Chapman  Model  of  Nighttime  Electron  Density  Profiles 
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Figure  2-3.  Contours  of  Zero-Degree  Aspect  Angle  for  10  MHz  E-Layer  Echoes 
as  Viewed  from  a Midlatitude  Location,  with  Feldstein-Starkov 
Auroral  Belt  Model  for  Various  Levels  of  Magnetic  Activity  at 
1200  Hours  Local  Time,  Daytime  Ionospheric  Model  A 
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Figure  2-4,  Contours  of  \ 
as  Viewed  fri 
Auroral  Belt 
1200  Hours  I 
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Figure  2-4. 


Contours  of  Zero-Degree  Aspect  Angle  for  20  MHz  E-Layer  Eehoes 
as  Viewed  from  a Midlatitude  Location,  with  Feldstein-Starkov 
Auroral  Belt  Model  for  Various  Levels  of  Magnetic  Activity  at 
1200  Hours  Local  Time,  Daytime  Ionospheric  Model  A 
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Figure  2-5.  Contours  of  3 
as  Viewed  fn 
Auroral  Belt 
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Contours  of  Zero-Degree  Aspect  Angle  for  30  MHz  E- Layer  Echoes 
as  Viewed  from  a Midlatitude  Location,  with  the  Feldsteln-Starkov 
Auroral  Belt  Model  for  Various  Levels  of  Magnetic  Activity  at 
1200  Hours  Local  Time,  Daytime  Ionospheric  Model  A 
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Figure  2-6. 
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gure  2-6. 


Contours  of  Zero-Degree  Aspect  Angle  for  10  MHz  E- Layer  Echoes 
as  Viewed  from  a Midlatitude  location,  with  the  Feldstein-Starkov 
Auroral  Belt  Model  for  Various  Levels  of  Magnetic  Activity  at 
1200  Hours  Local  Time,  Daytime  Ionospheric  Model  B 
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Figure  2-7.  Contours  of  Ze 
as  Viewed  fron 
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Figure  2-7,  Contours  of  Zero-Degree  Aspect  Angle  for  20  MHz  E-Layer  Echoes 
as  Viewed  from  a Midlatitude  Location,  with  the  Feldstein-Starkov 


Auroral  Belt  Model  for  Various  Levels  of  Magnetic  Activity  at 


1200  Hours  Local  Time,  Daytime  Ionospheric  Model  B 
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Figure  2-8.  Contours  of  Zero- 
as  Viewed  from  a 
Auroral  Belt  Modi 
1200  Hours  Local 
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Figure  2-8.  Contours  of  Zero-Degree  Aspect  Angle  for  30  MHz  E-Layer  f^choes 
as  Viewed  from  a Midlatitude  Location,  with  the  Feldstein-Starkov 
Auroral  Belt  Model  for  Various  Levels  of  Magnetic  Activity  at 
1200  Hours  Local  Time,  Daytime  Ionospheric  Model  B 


It  is  evident  that  the  spatial  distribution  of  the  perpendicular  contour  lines  is  a 
function  of  the  transmission  frequency  and  the  electron  density  profile.  In  addition,  it  is 
seen  that  it  is  also  a function  of  the  azimuth-elevation  angle  orientation  which,  in  turn, 
specifies  the  magnetic  inclination  and  declination  angles  in  space.  Magnetic  field  ortho- 
gonality can  be  attained  on  both  the  upward  and  downward  ray,  i.  e.  , after  the  ray  has  under- 
gone ionospheric  reflection.  In  fact,  normality  can  be  reached  at  multiple  altitudes 
(Millman,  1975). 

As  a first  approximation,  it  is  assumed  that  radar-auroral  echoes  can  be  expected  to 
appear  in  the  areas  formed  by  the  intersection  of  the  Feldstein-Starkov  auroral  ovals  and 
the  spatial  positions  where  the  magnetic  field  orthogonality  is  attained.  The  geographic 
coordinates  of  the  boundaries  of  the  Feldstein-Starkov  auroral  ovals  were  deduced  by  the 
nomograph  method  developed  by  Whalen  (1970). 

It  is  apparent  that  daytime  E-layer  HF  radar -auroral  echoes  should  not  be  observed. 
However,  as  shown  in  Figures  2-9  and  2-10,  E-layer  auroral  echoes  could  be  present  during 
the  nighttime  up  to  approximately  780  nmi  range  in  a northwesterly  direction. 

2.4  F-LAYER  BACKSCATTER  ECHO  REGIONS 

Contour  plots  of  where  spread-F  was  present  40,  60  and  80  percent  of  the  time  at 
1300  hours  local  standard  time  for  July  and  September  1957  (Penndorf,  1962)  are  shown  in 
Figures  2-11  through  2-16  together  with  magnetic  field  perpendicular  contour  lines.  The 
area  formed  by  the  intersections  of  the  spread-F  regions  and  the  normality  contours  defines 
the  region  where  spread-F  echoes  are  expected  to  occur. 

It  is  evident  that  the  40  percent  contour  lines  do  not  intersect  the  orthogonality  regions 
which  would  indicate  that  the  probability  of  spread-F  echoes  existing  during  the  daytime 
should  be  less  than  40  percent. 

The  occurrence  of  radar  echoes  from  the  F-layer  can  also  be  examined  in  terms  of 
the  irregularity  region  which  was  derived  from  scintillation  measurements  by  Aarons 
(1973b).  As  shown  in  Figures  2-17  and  2-18,  echoes  from  F-layer  irregularities  during 
the  daytime  could  exist  over  a wide  azimuthal  extent  and  to  ranges  out  to  over  1200  nmi. 


2-19/2-20 


UPWARD  RAY 
■DOWNWARD  RAY 
E • ELEVATION  ANGLE 


Figure  2-9.  Contours  of  Zero- 
as  Viewed  from  a 
Auroral  Belt  Mod* 
0000  Hours  Local 
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'igure  2-9.  Contours  of  Zero-Degree  Aspect  Angle  for  10  MHz  L-  y 

as  Viewed  from  a Midlatitude  Location,  with  the  Feldstein-Starkov 
Auroral  Belt  Model  for  Various  Levels  of  Magnetic  Activity  at 
0000  Hours  Local  Time,  Nighttime  Ionospheric  Model  A 
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Figure  2-10.  Contoui 
as  Viev 
Aurora: 
0000  He 
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Figure  2-10.  Contours  of  Zero-Degree  Aspect  Angle  for  10  MHz  E-Layer  Echoes 
as  Viewed  from  a Midlatitude  Location,  with  the  Feldstein-Starkov 
Auroral  Belt  Model  for  Various  Levels  of  Magnetic  Activity  at 
0000  Hours  Local  Time,  Nighttime  Ionospheric  Model  B 
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Figure  2-11.  Contours  of  Zero-Degree  Aspect  Angle  at  10  MHz  and  Penndorf's 
Probability  of  Occurrence  of  Spread-F  as  Viewed  from  a Mid- 
latitude Location,  at  1300  Hours  Local  Time,  Daytime  Ionospheric 
Model  A 
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Figure  2-12.  Contours  of  Zero-Degree  Aspect  Angle  at  20  MHz  and  Penndorf's 
Probability  of  Occurrence  of  Spread-F  as  Viewed  from  a Mid- 
latitude Location,  at  1300  Hours  Local  Time,  Daytime  Ionospheric 
Model  A 
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Figure  2-13.  Contours  of  Zero-Degree  Aspect  Angle  at  30  MHz  and  Penndorf’s 
Probability  of  Occurrence  of  Spread-F  as  Viewed  from  a Mid- 
latitude Location,  at  1300  Hours  Local  Time,  Daytime  Ionospheric 
Model  A 


Figure  2-14.  Contours  of  Zero-Degree  Aspect  Angle  at  10  MHz  and  Penndorf's 
Probability  of  Occurrence  of  Spread-F  as  Viewed  from  a Mid- 
latitude Location,  at  1300  Hours  Local  Time,  Daytime  Ionospheric 
Model  B 
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Figure  2-15.  Contours  of  Zero-Degree  Aspect  Angle  at  20  MHz  and  Penndorf's 
Probability  of  Occurrence  of  Spread-F  as  Viewed  from  a Mid- 
latitude Location,  at  1300  Hours  Local  Time,  Daytime  Ionospheric 
Model  B 
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Figure  2-16.  Contours  of  Zero-Degree  Aspect  Angle  at  30  MHz  and  Penndorf's 
Probability  of  Occurrence  of  Spread-F  as  Viewed  from  a Mid- 
latitude Location,  at  1300  Hours  Local  Time,  Daytime  Ionospheric 
Model  B 
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Figure  2-17.  Contours  of  Z 
Irregularity  R 
at  1200  Hours 


Figure  2-17.  Contours  of  Zero  Degree  Aspect  Angle  and  Aarons’  F-Layer 
Irregularity  Region  as  Viewed  from  a Midlatitude  Location, 
at  1200  Hours  Local  Time,  Daytime  Ionospheric  Model  A 
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Figure  2-18.  Contours  of  Zero-Degree  Aspect  Angle  and  Aarons’  F- Layer 
Irregularity  Region  as  Viewed  from  a Midlatitude  Location, 
at  1200  Hours  Local  Time,  Daytime  Ionospheric  Model  B 
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Fiffu  res  2-19  and  2-20  are  plots  of  the  zero  degree  aspect  angle  at  10  MHz  and 
Penndorf's  (19G2)  spread-F  probability  contours  at  0100  hour  local  time  for  both  nighttime 
models.  It  is  seen  that  the  magnetic  field  orthogonality  lines  intersect  the  HO  percent 
probability  contours.  This  would  imply  that  echoes  from  field-aligned  spread-F  should  be 
readily  observed  during  the  nighttime. 

Figure  2-21  which  depicts  Aarons'  (1973b)  F-layer  nighttime  irregularity  region  with 
the  zero  degree  aspect  angle  contours  reveals  that,  even  during  low  magnetic  activity 
Kp  = 0,  1),  an  II F radar  should  detect  reflections  from  F-layer  irregularities. 

2.  5 SUMMARY 

Tables  2-2  and  2-3  summarize  the  regions  (azimuth  extent  - A A,  and  range  extent  - 
A R)  where  field-aligned  HF  echoes  are  predicted  for  a midlatitude  location. 

It  is  seen  that  daytime  E -layer  auroral  echoes  should  be  nonexistent  while,  during 
the  nighttime,  for  relatively  quiet  magnetic  conditions  (Q  = 3),  the  echoes  could  appear 
within  an  azimuthal  spread  between  -55^  to  +20°  and  range  extent  from  280  to  780  nmi. 

Utilizing  Penndorf's  spread-F  probability  contours,  spread-F  echoes  are  predicted 
to  be  present  only  during  the  nighttime  at  the  low  frequencies  (in  the  vieinity  of  10  MHz). 

In  terms  of  Aarons'  irregularity  region,  however,  both  daytime  and  nighttime  F-layer 
echoes  can  be  expected. 

It  is  emphasized  that  the  regions  in  space  where  orthogonality  with  the  magnetic 
field  is  attained  is  a function  of  the  electron  density  along  the  ray  paths.  Thus,  the 
estimates  of  the  range  and  azimuth  orientation  of  the  field-aligned  echoes,  tabulated  in 
Tables  2-2  and  2-3,  will  be  modified  depending  on  the  ionization  distribution  in  the 
ionosphere. 

HF  baekseatter  radar  observations  conducted  at  Caribou,  Maine  (geographic  coordin- 
ates: 47°N,  68°W)  revealed  that  radar  reflections  from  E-layer  irregularities  could  at  times 
be  observed  over  an  entire  90°  azimuthal  coverage.  The  range  extent  of  E-layer  echoes 
attained  a value  as  large  as  500  km  (Campbell  et  al. , 1972).  The  angular  extent  of  the 
baekseatter  from  F-layer  irregularities  varied  from  1°  to  about  30°  while  the  range  extent 
was  found  to  be  as  great  at  2000  km.  The  height  interval  over  which  the  echoes  existed 
could  not  be  determined  because  of  equipment  limitations. 
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Figure  2-19.  Contours  of  Zero-Degree  Aspect  Angle  at  10  MHz  and  Penndorf's 
Probability  of  Occurrence  of  Spread-F  as  Viewed  from  a Mid- 
latitude Location,  at  0100  Hours  Local  Time,  Nighttime  Ionospheric 
Model  A 
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Figure  2-20.  Contours  of  Zero-Degree  Aspect  Angle  at  10  MHz  and  Penndorf's 
Probability  of  Occurrence  of  Spread-F  as  Viewed  from  a Mid- 
latitude Location,  at  0100  Hours  Local  Time,  Nighttime  Ionospheric 
Model  B 
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Figure  2-21. 


Contours  of  Zero-Degree  Aspect  Angle  and  Aarons’  F-Layer 
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2 G MAGNETIC  ACTIVITY  INDICES 


It  should  bo  noted  that  the  magnetic  activity  index,  Q,  used  by  I’eldstein  and  Starkov 
Q9G7)  is  related  to  the  geomagnetic  planetary  index,  Kp.  According  to  Valley  (19G5),  a 
Q 0 is  comparable  to  a Kp  = 0 or  1;  Q 1,  Kp  = 2;  Q - 4,  Kp  - 5 or  G and  Q 7, 

Kp  = 8 or  9. 

The  geomagnetic  planetary  index,  Kp,  is  a mean  three-hour  magnetic  reading  derived 
from  observations  located  between  47°  and  G3°  geomagnetic  latitude.  The  scale  of  0 to  9 is 
an  arbitrary  one  where  0 refers  to  very  quiet  magnetic  conditions  End  9 to  extremely 
disturbed  magnetic  conditions. 

The  magnetic  activity  index,  Q (often  referred  to  as  the  Polar  Range  Index),  is  a 
quarter -hourly  index  which  is  at  times  used  to  indicate  the  severity  of  magnetic  conditions 
for  high  latitudes  (above  about  58°  geomagnetic  latitude).  The  index  has  not  been  widely 
employed  in  geophysical  analysis  because  of  the  complexity  involved  in  the  data  reduction 
of  magnetic  observations  to  obtain  the  parameter. 

The  cumulative  distribution  function  of  the  occurrence  of  Kp  is  plotted  in  Figure  2-22 
for  the  two  extremes  of  solar  activity  and  for  the  1932-1970  period  (Zawalick  and  Cage, 
1971).  In  the  period  from  April  1957  through  December  1958,  which  encompasses  the 
International  Geophysical  Year  (IGY),  the  sunspot  number  was  approximately  200  while, 
from  January  1964  through  October  1965,  which  was  a period  of  quiet  solar  activity,  the 
sunspot  number  was  10.  It  is  seen  that,  during  a time  of  minimum  solar  activity,  the 
median  and  upper  decile  value  of  Kp  is  1 and  3,  respectively,  while,  during  maximum  solar 
activity,  the  corresponding  values  increase  to  approximately  2 and  4,  respectively. 

Figure  2-23  is  a plot  of  the  occurrence  of  the  geomagnetic  planetary  index,  Kp,  from 
1950  through  1971  (Aarons,  1973a)  and  the  mean  sunspot  number.  It  is  apparent  that  the 
occurrence  of  magnetic  disturbed  conditions  is  not  correlated  with  solar  activity.  These 
results  somewhat  contradict  the  statistical  data  presented  in  Figure  2-22. 
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SECTION  III 


SIGNAL  AMPLITUDE  AND  CROSS  SECTION  STATISTICS  OF 
IONOSPHERIC  BACKSCATTER  ECHOES 


3.  1 INTRODUCTION 

In  this  section,  an  analysis  is  made  of  the  statistical  characteristics  of  the  amplitude 
and  cross-sectional  area  of  HF  ionospheric  backscatter  echoes  that  could  be  observed 
by  an  HF  backscatter  radar  operating  in  the  midlatitudes  with  a northerly  oriented  antenna 
beam.  The  echo  data  that  are  evaluated  were  obtained  from  measurements  made  by  the 
Stanford  Research  Institude  (SRI)  (Leadabrand  et  al.f  1965)  and  under  the  Polar  Fox  II 
program  (Bradley  et  al. , 1972;  Herman  and  Vargas-Vila,  1973). 

3.  2 RADAR  EQUATION 

In  order  to  interpret  the  SRI  and  Polar  Fox  II  auroral  data  in  terms  of  the  correspond- 
ing effects  on  an  HF  backscatter  radar,  it  is  necessary  to  determine  the  relative  sensitivities 
of  the  radars  to  the  aurora.  This  is  accomplished  by  comparing  the  various  parameters 
which  enter  into  the  radar  equation. 


According  to  the  radar  equation,  the  signal -to -noise  ratio  (SNR)  of  the  received  signal 
is  given  by 


P 

r 
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P.G.G  A <7 
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(4tt)3  R4  KTB  NF  L L 
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where  Pt  is  the  transmitted  power,  Gt  is  the  gain  of  the  transmitting  antenna,  Gr  is  the  gain 

of  the  receiving  antenna,  A is  the  wavelength,  <r  is  the  radar  cross  section  of  the  target, 
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R is  the  radar  range,  K is  Boltzmann’s  constant  (1.  38  x 10  Ws/°K),  T is  the  ambient 

temperature  (288°K),  B is  the  receiver  noise  bandwidth,  NF  is  the  receiver  noise  figure, 

L is  the  system  loss  and  L is  the  two-way  loss  due  to  the  propagation  medium.  This  re- 
s p 

lationship  is  applicable  to  a point  target,  i.  e. , the  antenna  beamwidth  is  much  larger  than 
the  dimensions  of  the  target. 
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For  aurora  which  partially  fills  the  antenna  beam,  <j  is  replaced  by  the  volume 
scattering  coefficient,  c r , i.  o, , radar  cross  section  per  unit  volume,  such  that 

V 

a = <tv  V (3-2) 

where  V is  the  volume  filled  by  the  scatterers. 

The  volume  of  the  aurora  Intersected  by  the  antenna  beam,  shown  in  Figure  3-1,  can 
be  represented  by 

V = R ft  AR  Ah  (3-3) 

where  ft  is  the  horizontal  beamwidth  of  the  antenna,  Ah  is  the  thickness  of  the  aurora  in 
the  vertical  direction.  The  parameter,  AR  is  the  radar  range  resolution  defined  by 

AR  ~ (3-4) 


where  c is  the  free-spacc  velocity  and  t is  the  pulse  length. 

Substituting  Equations  (3-2),  (3-3)  and  (3-4)  in  Equation  (3-1),  it  follows  that 

P P,  G,  G A2  cr  cr  ft„  Ah 
r _ t t r v H 

N 2(4tt)3  R3  KTB  NF  L L 
v ' s p 

It  is  seen  that,  for  a partially  filled  beam,  the  SNR  is  inversely  proportional  to  the 
range  cubed. 

When  the  scatterers  completely  fill  the  antenna  beam,  the  scattering  volume  for  a 
rectangular  antenna  can  be  described  by 


(3-5) 


V = R Pu  Pv  AR 

where  fty  is  the  vertical  antenna  beamwidth. 

Substituting  Equations  (3-2),  (3-4)  and  (3-6)  in  Equation  (3-1)  results  in 


(3-6) 


^ = Pt  Gt  Gr  ^ % CT 
N 2(4tt)3  R2  KTB  NF  Lg  Lp 


0-7) 


It  is  of  interest  to  note  that,  for  the  filled  beam  case,  the  SNR  is  inversely  proportional 
to  the  range  squared. 
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BEST  AVAILABLE  COPY 


Figure  3-1.  Intersection  of  Antenna  Beam  with  Auroral  Arc 
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3.  3 SHI  E-LAYEH  AURORAL  DATA 


During  the  1959-GO  period,  SRI  conducted  an  experimental  program  at  Fraserburgh, 
Scotland,  for  the  study  of  the  characteristics  of  radar-auroral  echoes  (Leadabrand,  et  al. , 
1965).  Radar  measurements  were  made  simultaneously  at  a frequency  of  30,  401  and  800 
MHz.  The  amplitude  of  the  auroral  echoes  that  could  be  encountered  by  a high-powered  HF 
baekscatter  radar  is  obtained  by  extrapolation  from  the  SRI  data. 

Figure  3-2  is  a plot  of  the  relative  signal  strengths  of  auroral  echoes  recorded  at  the 
three  frequencies  in  April  1960.  The  region  of  uncertainty  between  30  and  401  MHz  results 
from  the  fact  that  there  was  no  way  to  specify  to  what  degree  the  auroral  seatterers  filled 
the  30-MHz  antenna  beam.  At  401  and  800  MHz,  the  antenna  beamwidth  was  1.  2°  while,  at 
30  MHz,  it  was  15°.  Thus,  it  was  assumed  that  the  seatterers  completely  filled  the  beam 
at  the  two  higher  frequencies.  However,  it  was  doubtful  that  the  beam  was  completely 
filled  at  30  MHz.  The  1°  x 1°  aurora  can  be  considered  to  be  a point  target  in  terms  of  the 
30-MHz  antenna. 

It  is  noted  that  the  data  in  Figure  3-2  are  normalized  with  respect  to  401  MHz.  For  a 
filled  beam  antenna,  the  relative  signal  strength  of  the  30-MHz  echoes  was  6 dB  greater  than 
those  at  401  MHz,  while  for  the  point  target  case,  there  was  a 28-dB  difference  between 
the  30-  and  401-MHz  data. 

The  absolute  magnitude  of  the  30-MHz  auroral  echoes  can  be  deduced  by  comparing 
the  system  sensitivities  of  the  30-  and  401-MHz  radars.  According  to  Table  3-1,  the 
401-MHz  radar  was  48.  9-dB  more  sensitive  than  the  30-MHz  radar  assuming  that  the 
aurora  was  a point  target,  and  26.  9-dB  more  sensitive,  assuming  that  the  30-MHz  antenna 
beam  was  completely  filled  with  seatterers.  Thus,  as  indicated  in  Table  3-2,  it  follows  that 
the  magnitude  of  the  30-MHz  radar  auroral  data  was  20.  9 dB  less  than  that  of  the  401-MHz 
data. 

The  amplitude  distribution  of  the  401-MHz  radar -auroral  echoes  is  presented  in 
histogram  form  in  Figure  3-3.  It  is  noted  that  the  system  noise  level  was  approximately 
-133  dBm  and  that  the  system  saturated  at  -60  dBm  which  accounts  for  the  large  number 
of  echoes  at  that  amplitude. 
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Figure  3-2.  Relative  Signal  Strength  of  Auroral  Echoes  at  30-,  401-  and  800-MHz 

Obtained  by  the  SRI  Radars  at  Scotland  (After  Leadabrand,  et  al.  , 1965) 


TABLE  3-1 


COMPARISON  OF  SRI-SCOTLAND  30-MIIz  AND  401-MIIz  RADAR  SYSTEM 
SENSITIVITIES  UTILIZED  IN  THE  APRIL  19G0  RADAR-AURORAL  OBSERVATIONS 


Parameter 

30 -MHz  Radar 

401-MHz  Radar 

Point  Target 
(401  MHz/30  MHz) 
(dB) 

Scatterers  Fill 
Antenna  Beam 
(401  MHz/30  MHz) 
(dB) 

2 2 
A (m  ) 

(10.  o)2 

(0.  748) 2 

-22.  5 

-22.  5 

P (MW) 

0.  0015 

0.  12 

+19 

+19 

Gt  Gr  (dB) 

40 

84 

+44 

+44 

(BW)a  (deg) 

15 

1.  2 

- 

-11 

(BW)e  (deg) 

15 

1.2 

- 

-11 

T (US) 

300 

300 

0 

0 

B (kHz) 

6 

G 

0 

0 

NF  (dB) 

9 

4.5 

+4.5 

+4.5 

Lg  (dB) 

7 

3.  1 

+3.  9 

+3.  9 

+48.  9 

+26.9 

Note:  (BW)a  = Azimuth  Beamwidth 

<BW)  = Elevation  Beamwidth 
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TABLE  3-2 

ESTIMATE  OF  SRI  30-MlIz  RADAR-AURORA L SIGNA L MAGNITUDE  RELATIVE  TO 

401 -MHz  RADAR  DATA 


Normalized 
Auroral  Signal 

Radar 

Auroral  Signal 

Auroral  Scatterers- 

M gnitude 

Sensitivities 

Magnitude 

Antenna  Beam 

(30  Mriz/401  MHz) 

(401  MHz/30  MHz) 

(30  MHz/401  MHz) 

Configuration 

(dB) 

(dB) 

(dB) 

Point  Target 

+28 

+48.9 

-20.  9 

Filled  Beam 

+6 

+26.9 

-20.  9 

SATURATION 


Figure  3-3.  Amplitude  Histogram  of  SRI  401-MHz  Scotland  Radar  Auroral  Data 
(After  Leadabrand,  et  al. , 1965) 
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The  histogram  data,  when  replotted  as  a cumulative  distribution  on  probability 
paper,  shown  in  Figure  3-4,  closely  fit  a straight  line.  This  characteristic  is  indicative 
that  the  401-MHz  amplitudes  are  log-normally  distributed.  The  theoretical  normal  distri- 
bution function  plotted  in  Figure  3-3  was  computed  utilizing  the  value  of  -77.  2 dBm  for  the 
mean  as  obtained  from  Figure  3-4. 

It  is  assumed  in  this  analysis  that  the  30-MHz  amplitude  data  also  followed  a 
Gaussian  distribution.  Table  3-3  listis  the  statistical  parameters  of  the  401-MHz  distribu- 
tion and  the  extrapolated  30-MHz  distribution.  The  latter,  which  is  also  shown  in  Figure 
3-4,  was  derived  on  the  basis  that  the  30-MHz  auroral  echoes  were  20.  9 dB  less  than  the 
401-MHz  results. 

In  estimating  the  magnitude  of  the  auroral  echoes  that  could  be  observed  by  an  HF 
backscatter  radar,  it  is  necessary  to  compare  the  sensitivity  of  the  HF  baekscattcr  radar 
to  that  of  the  SRI  30-MHz  radar.  According  to  Table  3-4,  assuming  a pulse  length  of 
400  ps,  the  HF  radar  is  37. 1 dB  more  sensitive  than  the  SRI  30-MHz  radar  for  the  point 
target  case  and  33.  5 dB  more  sensitive  for  the  filled  beam  case.  When  the  pulse  length  is 
decreased  to  10  ps,  the  difference  in  sensitivities  decreases  to  21. 1 dB  and  1.  5 dB  for  the 
point  target  and  filled  beam  case,  respectively. 

The  predicted  HF  backscatter  radar -auroral  clutter  levels  at  30  MHz  are  presented 
in  Table  3-5.  The  calculations  are  based  on  combining  the  SRI  30  MHz  radar-auroral  data 
given  in  Table  3-3  with  radar  sensitivities  in  Table  3-4. 

In  deducing  the  radar-auroral  signal  levels  at  other  frequencies  in  the  HF  band,  it 
is  assumed  that  the  auroral  echo  power  is  frequency  dependent  according  to  the  law 

P = kf"n  (3-8) 

where  k and  n are  constants.  Utilizing  the  data  in  Figure  3-2,  it  can  be  shown  that,  for 
this  power  law,  n = 2.  5 for  a point  target  and  n = 0.  5 for  the  filled  beam  case.  However, 
an  examination  of  other  radar-auroral  data  taken  by  SRI  at  Scotland  indicated  a value  of 
n = 6.  8 for  a point  target  and  n = 2.  7 for  a filled  beam  (Leadabrand,  et  al. , 1965). 

An  analysis  of  auroral  echoes  from  simultaneous  multiple  frequency  observations 
in  Alaska  by  Leadabrand  et  al.  (1967)  revealed  that  n = 2 for  frequencies  between  50  and 
850  MHz  and  n = 5 between  850  and  3000  MHz. 
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Figure  3-4.  Cumulative  Distributi 
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B4509 


Function  of  SRI  Scotland  Radar -Auroral  Data 


TABLE  3-3 

STATISTICAL  PARAMETERS  OF  SRI  RADAR-AURORA  L DATA 


Statistical 

Parameter 

401-MHz  Data 
(dBm) 

30-MHz  Data 
(dBm) 

Upper  Decile 

-58.  2 

-79.1 

Upper  Quartile 

-67.  2 

1 

00 

00 

h-» 

Median 

-77.  2 

-98.1 

Lower  Quartile 

-87.  2 

-108.1 

Lower  Decile 

-96.2 

-117. 1 

Standard  Deviation 

14.8 

14.8 

TABLE  3-4 

COMPARISON  OF  IIF  BACKSCATTER  RADAR  SYSTEM  SENSITIVITY 
1 O THAT  OF  THE  SRI-SCOTLAND  30-MIIz  RADAR 


Parameter 


P (MW) 


IIF  Backs eatter 
Radar 

SRI-Scotland 

Radar 

Point  Target 
(IIF  Radar /SRI) 
(0B) 

Scatterers  Fill 
Antenna  Beam 
(HF  Radar /SRI) 
(dB) 

0 

9 

(10.  Of 

(10.  Of 

0 

0 

0.  4 

0.  0015 

+24.3 

+ 1.3 

45 

40 

+ 5 

+5 

5 

15 

- 

-4.8 

15 

15 

- 

0 

400 

300 

_ 

+1.2 

10 

-14.  8 

2.5 

6 

+3.  8 

+3.  8 

100 

-12.  2 

-12.  2 

9 

9 

0 

0 

3 

7 

+4 

+4 

(T  = 400  ps) 

+37.  1 dB 

+33.  5 dB 

(t  - io  ns) 

1 

+21.  1 dB 

+ 1.  5 dB 

I 


TABLE  3-5 

ESTIMATE  OF  HF  BACKSCATTER  RADAR-AURORAL  SIGNAL  LEVELS 
AT  30  MHz  BASED  ON  SRI  DATA 


Point  Target 

Scatterers  Fill 
Antenna  Beam 

Statistical 

Parameter 

t = 400  ps 
(dBm) 

T = 10  pS 
(dBm) 

t = 400  ps 
(dBm) 

T = 10  PS 
(dBm) 

Upper  Decile 

-42.  0 

-58.  0 

-45.6 

-77.6 

Upper  Quartile 

-51.0 

-67.0 

-54.6 

-86.6 

Median 

-61.  0 

-77.  0 

-64.6 

-96.6 

Lower  Quartile 

-71.  0 

-87.  0 

-74.6 

-106.6 

Lower  Decile 

-80.  0 

-96.  0 

-83.6 

-115.  6 

Standard  Deviation 

! 

14.8 

14.  8 

14.8 

14.8 

Radar  backscattor  measurements  of  artificial  electron  clouds  in  the  E-region  of  the 
ionosphere  by  Gallagher  and  Barnes  (19G3)  yielded  a constant,  n,  of  4 for  frequencies  be- 
tween approximately  20  and  50  MHz  and  -4  between  5 and  20  MHz.  At  times,  however,  it 
was  found  that  the  amplitude  returns  were  insensitive  to  frequency  in  approximately  the 
10-  to  20-Mllz  range. 

From  simultaneous  auroral  echo  measurements  at  a frequency  of  40.7,  143.5  and 
22G  MHz,  Flood  (19G5)  has  deduced  a value  of  n = 3.  5 to  3.  7 between  49.  7 and  143.  5 MHz 
and  n = G.  5 ± 1.  3 between  143.  5 and  22G  MHz.  During  periods  of  intense  auroral  reflections, 
there  were  instances  when  n = 12.  9 ± 1.  2 between  143.  5 and  22G  MHz. 

Since  there  is  wide  discrepancy  in  the  experimental  measurements  of  the  frequency 
dependence  of  auroral  backscatter,  which  could  be  due  to  the  characteristics  of  the  auroral 
ionization,  i.  e.  , inhomogeneous  distribution  of  auroral  electrons,  varying  scale  sizes  of 
ionization  irregularities  and  different  scattering  altitudes,  a value  of  4 for  the  exponent  in 
Equation  (3-8)  is  assumed  in  this  analysis. 

The  estimated  median  auroral-eeho  amplitudes  that  could  confront  an  HF  baekseatter 
radar  in  the  5-  to  30-MHz  band  are  plotted  in  Figure  3-5.  The  external  noise  level  of  a 
rural  environment  also  shown  in  Figure  3-5  is  the  average  median  value  in  a 1-Hz  bandwidth 
as  predicted  for  a location  in  Maine  for  all  seasons,  all  times  of  day  and  a sunspot  number 
of  70.  The  more  precise  values  of  the  external  noise  levels  are  given  in  Table  3-6.  It 
is  seen  that  the  auroral  echo  amplitude  increases  with  decreasing  frequency  and  increasing 
pulse  length  and  that  the  amplitude  of  the  point  target  case  is  greater  than  that  of  the  filled 
beam  case. 

Figure  3-6  is  a replot  of  the  data  in  Figure  3-5  in  terms  of  the  auroral  SNR  assum- 
ing a 1-Hz  bandwidth.  It  is  of  interest  to  note  that  the  SNR  maximizes  in  the  vicinity  of 
10  MHz.  This  is  due  to  the  fact  that  the  external  noise  level  increases  with  decreasing 
frequency. 

The  statistical  estimates  of  the  auroral  SNR  are  listed  in  Table  3-7.  It  is  seen 
that  the  quartile  values  are  ±10  dB  with  respect  to  the  median  level  while  the  deeile  values 
are  ±19  dB.  In  other  words,  there  is  a 20-dB  difference  in  the  auroral  elutter-to-noise 
ratio  between  the  upper  and  lower  quartiles  and  38-dB  difference  between  the  upper  and 
lower  deciles. 
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Figure  3-5.  Estimated  Median  Radar -Auroral  Signal  Amplitude  for  an  HF 
Backscatter  Radar  Based  on  SRI  Scotland  Radar  Data 
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Figure  3-6.  Estimated  Median  Radar -Auroral  Signal-to-Noise  Ratio  for  an  HF 
Backscatter  Radar  Based  on  SRI  Scotland  Radar  Data 


TABLE  3-6 


AVERAGE  MEDIAN  VALUES  OF  EXTERNAL  NOISE  (RURAL  ENVIRONMENT) 
AS  A FUNCTION  OF  MEDIAN  MAXIMUM  USABLE  FREQUENCY  AS 
PREDICTED  FOR  MAINE  FOR  ALL  SEASONS, 

TIMES  OF  DAY  AND  SSN  OF  70 


Maximum  Usable 
Frequency 
(MHz) 

Median  Noise 
(dBm/Hz) 

5 

-122.6 

10 

-136.8 

15 

-142.  6 

20 

-146.5 

25 

-148.5 

30 

-150.  3 

TABLE  3-7 

ESTIMATE  OF  HF  BACKSCATTEK  RADAR  E-LAYER  AURORAL  SNR 

BASED  ON  SRI  DATA 


Pulse 

Length 

(us) 


Upper 

Decile 

(dB) 


111.  7 
95.  7 
108.  1 

76.  1 

113.  9 
97.  9 

110.  3 
78.  3 

112.6 

96.6 
109.  0 

77.  0 

111.  5 
95.  5 

107.  9 
75.  9 

109.  7 

93.7 
106.  1 

74.  1 

108.  3 
92.  3 

104.  7 

72.7 


Upper 

Quartile 

(dB) 


103.  6 
87.6 
100.  0 
68.  0 


102.5 
86.  5 
98.  9 
66.  9 

100.  7 
84.  7 
97.  1 
65.  1 


Median 

(dB) 


Lower 

Quartile 

MB) 


Lower 

Decile 

(dB) 


PT  = Point  Target 
FB  = Filled  Beam 


Standard  Deviation  = 14.  8 dB 


The  auroral  SNR's  presented  in  Table  3-8  are  the  averages  of  the  filled  beam  and  the 
point  target  echo  data  of  Table  3-7.  The  interesting  point  to  note  in  Table  3-8  is  that  the 
averaging  manipulation  does  not  alter  the  value  of  the  differences  between  the  magnitudes 
of  the  statistical  parameters. 

It  should  be  noted  that  the  height  distribution  of  the  auroral  echoes  observed  by  SRI  in 
Scotland  were  peaked  at  about  100  to  120  km,  i.  e. , in  the  E-region,  although  heights  as 
great  as  200  km  were  observed.  Thus,  the  estimates  given  in  Tables  3-7  and  3-8  can  be 
considered  to  apply  only  to  E -region  reflections. 

An  additional  item  which  should  be  mentioned  is  that  the  amplitude  distribution  of  the 
401-MHz  auroral  echoes,  shown  in  Figure  3-3,  which  was  used  as  a basis  of  extrapolation 
to  the  HF  backscatter  radar,  was  applicable  to  data  obtained  when  the  direction  of  propaga- 
tion at  the  reflection  point  was  perpendicular  to  the  magnetic  field  lines.  When  the  angle 
between  the  directions  of  propagation  and  magnetic  field  lines  are  not  orthogonal,  the  auroral 
clutter  amplitude  should  decrease. 

Bates  and  Albee  (1969)  report  that,  if  a simple  model  of  auroral  backscatter  from  the 
E -region  that  includes  ionospheric  refraction  effects  is  assumed,  a lower  limit  of  6 dB/deg 
in  the  15-  to  50-MHz  range  was  observed  for  the  aspect-sensitive  decrease  in  backscatter 
cross  section  with  off-perpendicular  angle  from  the  geomagnetic  field. 

The  aspect-sensitive  decrease  of  the  cross  section  of  F-layer  backscatter  echoes  was 
found  by  Bates  and  Albee  (1970)  to  be  on  the  order  of  5 dB/deg  of  off-perpendicularity  from 
the  magnetic  field  in  the  6-  to  15-MHz  band. 

The  results  of  aspect  sensitivity  measurements  made  in  the  frequency  range  from  HF 
to  UHF  are  summarized  in  Table  3-9. 

3.4  POLAR  FOX  II  E -LAYER  AURORAL  AND  F-LAYER  DATA 

An  objective  of  the  Polar  Fox  II  program  was  to  investigate  the  characteristics  of  HF 
ionospheric  backscatter  echoes.  The  experimental  measurements  were  conducted  at 
Caribou,  Maine,  during  the  period  between  December  1971  and  November  1972.  Ionospheric 
observations  were  made  at  five  frequencies  in  the  HF  band.  However,  in  this  analysis,  only 
the  echo  data  recorded  at  8. 125  and  14.  875  MHz  in  January  1972  (Bradley,  et  al. , 1972)  and 
in  July  1972  (Herman  and  Vargas-Vila,  1973)  are  examined. 
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TABLE  3-8 


ESTIMATE  OF  HF  BACKSCATTER  RADAR-AVERAGE  E-LAYER  AURORAL  SNR 

BASED  ON  SRI  DATA 


Frequency 

(MHz) 

Pulse 

Length 

(MS) 

Upper 

Decile 

(dB) 

Upper 

Quartile 

(dB) 

Median 

(dB) 

Lower 

Quartile 

(dB) 

Lower 

Decile 

(dB) 

5 

400 

BB 

102.  0 

92.  0 

82.  0 

73.  0 

10 

EH 

86.0 

76.  0 

66.  0 

57.  0 

10 

400 

113.  2 

104.  2 

94.  2 

84.2 

75.  2 

10 

97.  2 

88.  2 

78.2 

68.  2 

59.  2 

15 

400 

111.9 

102.  9 

92.  9 

82.  9 

73.  9 

10 

95.  9 

86.  9 

76.9 

66.  9 

57.  9 

20 

400 

110.8 

101.  8 

91.8 

81.  8 

72.8 

10 

94.  8 

85.  8 

75.8 

65.  8 

56.  8 

25 

400 

109.0 

100.0 

90.  0 

80.0 

71.  0 

10 

93.  0 

84.  0 

74.0 

64.  0 

55.  0 

30 

400 

107.6 

98.6 

88.6 

78.6 

69.  6 

10 

91.6 

82.6 

72.6 

62.6 

53.  6 

I I 


Frequency 

(MHz) 

Aspect 
Sensi  vity 
(dB/deg) 

Aspect  Angle 
Interval 
(deg) 

Location 

Investigator 

G-15 

5 

0 < AO  < 5 

College,  Alaska 

Bates  and  Albee,  1970 

15-50 

>6 

0 <A 0 <4 

College,  Alaska 

Bates  and  Albee,  1969 

30 

2 

0 < AO  <1 

Fraserburgh, 

Leadabrand  et  al.  , 1965 

401,  800 

10 

0 < AO  < 7 

Scotland 

48 

2 

4.5  £ A 0 £19.  5 

Baker  Lake  and 
Saskatoon,  Canada 

McDiarmid,  1972 

48.  5 

1 

1°  £A0  £19° 

Baker  Lake, 
Ottawa,  Resolute 
and  Saskatoon, 
Canada 

McNamara,  1972 

50,  139,  398, 
850,  1210, 
3000 

10 

0 < AO  <2 

Homer,  Alaska 

Chesnut  et  al.  , 1968 

448 

5 

4.5  £A0£7.  5 

Prince  Albert, 
Canada 

Jaye  et  al.  , 1967 

501 

10 

0 <A0  <2 

Manchester, 

England 

Barber  et  al.  , 1962 

I 


Figures  3-7  and  3-8  are  plots  of  the  diurnal  variation  of  the  area  scattering  co- 
efficient of  E -layer  and  F-laycr  echoes,  respectively.  '1’he  latter  ean  he  considered  to 
originate  at  an  altitude  of  approximately  300  km  while  the  former  at  an  altitude  of  about 
100  km.  The  data  arc  averaged  over  the  azimuthal  extent  o;  the  radar,  i.  e.  , -30°  to  +00° 
measured  with  respect  to  geographic  north.  For  simplicity  in  data  analysis,  the  azimuthal 
spread  was  divided  into  5 zones  of  -30°  to  -12°  (zone  1),  -12°  to  +0°  (zone  2),  (-0°  to  +24° 
(zone  3),  +24°  to  +42°  (zone  4)  and  +42°  to  +00°  (zone  5). 


It  should  be  noted  that  the  data  in  Figures  3-7  and  3-8  describe  the  median  and  the 


lower  and  upper  quartile  values  and  that  the  propagation  losses,  L^,  have  not  been  removed 
from  the  calculation  of  the  area  scattering  coefficient,  <j..  In  other  words,  the  vertical 


axis  is,  in  reality,  the  term,  <r./L  . 

’ A p 


It  is  apparent  that  the  signals  from  both  E-  and  F-layer  irregularities  show  a slight 
decrease  in  the  median  value  of  the  scattering  coefficient  during  the  nighttime  hours.  How- 
ever, in  the  ease  of  the  8. 125-MIIz  F-layer  clutter,  the  median  scattering  coefficient  re- 


mains relatively  constant  at  approximately  -30  dBsm/m"  during  the  nighttime  and  decreases 

o 
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to  about  -40  dBsm/m  during  the  early  hours  of  the  morning.  In  addition,  it  is  seen  that,  at 
a given  time,  the  median  value  of  scattering  eoeficient  at  14.  875  MHz  is  less  than  that  at 
8.  125  MHz  irrespective  of  the  altitude  of  the  scattering  irregularities. 


The  zonal  variation  of  the  area  scattering  coefficient  of  the  E-  and  F-layer  irregular- 
ities are  presented  in  Tables  3-10  and  3-11,  respectively.  It  is  seen  that,  for  all  the  statis- 
tical parameters,  the  8. 125-MHz  data  are  from  12  to  21  dB  greater  than  the  14.  875-MHz 
results.  The  most  interesting  feature  of  the  data  is  the  lack  of  influence  of  the  geographic 
locations  of  the  seatterers  on  the  magnitude  of  the  scattering  eoeffieient.  By  intuition,  it 
wouldseem  that  the  scattering  coefficient  would  be  greater  in  zone  1 or  2 than  in  zone  4 or  5, 
sinee  the  former  are  the  regions  eneompassed  by  the  auroral  belt  (Feldstein  and  Starkov, 
1967). 


Table  3-12  illustrates  the  difference  between  the  F-layer  scattering  coefficient, 
(cr^)F,  and  the  E-layer  coefficient,  (cr^)^,  shown  in  Tables  3-10  and  3-11,  respectively. 
It  is  interesting  to  note  that,  except  for  8. 125-MHz  upper  quartile  results,  the  F-layer 
area  scattering  eoeffieient  is  on  the  order  of  17  to  20  dB  greater  than  the  E-layer 
eoeffieient,  irrespective  of  the  transmission  frequency. 
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LOCAL  TIME  (h) 


Figure  3-7.  Diurnal  Variation  of  Area  Scattering  Coefficient  of  Polar  Fox  II  E-Layer 
Irregularities  Averaged  over  Five  Zones  at  a Frequency  of  8. 125  and 
14.  875  MHz,  January  1972  (after  Bradley  et  al. , 1972) 


/L  (dBsm 


TABLE  3-10 

ZONAL  VARIATION  OF  AREA  SCATTERING  COEFFICIENT  OF  E-LAYER 
IRREGULARITIES  AVERAGED  OVER  TIME  AT  A FREQUENCY  OF 
8.  125  AND  14.  875  MHz,  JANUARY  1972 


Frequency 

(MHz) 

Statistical 

Parameter 

8.  125 

Upper 

Quartile 

Median 

Lower 

Quartile 

14.  875 

Upper 

Quartile 

Median 

IiOwer 

Quartile 

Area  Scattering  Coefficient  (dBsm/m  ) 


Zone 


2 

3 

-35.  8 

-37.  5 

-48.  2 

-51.  6 

-58.  8 

-62.  9 

-58.4 

-57.  0 

-66.2 

-66.  2 

-72.6 

-74.  2 

5 

Average 

-32.  5 

-34.  2 

-46.  0 

-47.6 

-56.  9 

-58.6 

-54.  9 

-56.5 

-62.  2 

-63.  9 

-74.3 

■ -74.0 

Note:  Propagation  losses  not  removed  from  area  scattering  coefficient. 


TABLE  3-11 


ZONAL  VARIATION  OF  AREA  SCATTERING  COEFFICIENT  OF  F-LAYER 
IRREGULARITIES  AVERAGED  OVER  TIME  AT  A FREQUENCY  OF 
8.  125  AND  14.  875  MHz,  JANUARY  1972 


2 

Area  Scattering  Coefficient  (dBsm/m  ) 

Frequency 

(MHz) 

Statistical 

Parameter 

Zone 

1 

2 

3 

4 

5 

Average 

8.  125 

Upper 

Quartile 

-28.  4 

-25.8 

-26.  8 

-28.3 

-30.  0 

-27.2 

Median 

-33.  5 

-32.2 

-33.  2 

-34.  8 

-36.6 

-33.  6 

Lower 

Quartile 

-42.  9 

-40.5 

-42.  1 

-45.  7 

-46.2 

-42.0 

14.875 

Upper 

Quartile 

-41.3 

-39.3 

-42.  8 

-40.  8 

-35.3 

-37. 1 

Median 

-49.  7 

-50.2 

-53.4 

-50.  1 

-45.4 

-47.2 

Lower 

Quartile 

-5G.  7 

-59.8 

-62.  2 

-63.  0 

-52.5 

-54.3 

Note:  Propagation  losses  not  removed  from  area  scattering  coefficient. 


TABLE  3-12 

DIFFERENCE  BETWEEN  AREA  SCATTERING  COEFFICIENTS  OF  F-LAYER  AND 
E-LAYER  IRREGULARITIES  AT  A FREQUENCY  OF  8.  125  MHz  AND 

14.  875  MHz,  JANUARY  1972 


Frequency 

(MHz) 

8. 125 


14. 875 


(°A^F  " ^A^E 
Zone 


Statistical 

Parameter 

1 

2 

3 

4 

5 

Average 

Upper 

Quartile 

8.5 

10.  0 

10.7 

9.9 

2.5 

9.6 

Median 

14.7 

16.  0 

18.4 

17.8 

9.4 

17.  1 

Lower 

Quartile 

1 

20.1 

18.3 

20.8 

17.2 

10.7 

19.5 

Upper 

Quartile 

17.9 

19. 1 

14.2 

18.3 

19.6 

18.5 

Median 

15.3 

16.0 

12.8 

18.0 

16.8 

16.6 

Lower 

Quartile 

18.7 

12.  8 

12.0 

13.6 

21.8 

20. 1 

Note:  Propagation  losses  not  removed  from  area  scattering  coefficients. 


The  conversion  of  the  area  scattering  coefficient  into  a radar  cross  section  can  bo 
accomplished  utilizing  Equations  (3-2),  (3-3)  and  (3-4).  It  is  evident  that  the  radar  cross 
section,  c r,  is  related  to  the  volume  scattering  coefficient,  a , by 

cr  = <rv  H ¥ Ah  (3_9) 

Since  the  thickness  of  the  scatterers  in  the  vertical  direction,  Ah,  was  an  unknown 
factor  in  the  Polar  Fox  II  experiment,  it  was  assumed  to  be  unity  (1  m).  Thus,  Equation 
(3-9)  can  be  described  in  terms  of  the  area  scattering  coefficient,  cr^,  by 

a = aA  A = aA  R /?H  ¥ (3"10) 

where  A is  the  area  containing  the  scatterers. 

The  scattering  area  of  the  HF  backscatter  radar  is  depicted  in  Figure  3-9  as  a 
function  of  range  for  the  two  pulse  lengths. 

The  cross  sections  of  ionospheric  backscatter  echoes  that  could  be  observed  by  the 
HF  backscatter  radar  are  presented  in  Table  3-13.  The  calculations  are  derived  utilizing 
the  data  in  Tables  3-4,  3-11  and  3-12  and  assuming  a radar  range  of  700  km  for  the  E- 
layer  irregularities  and  1500  km  for  the  F-layer  irregularities.  These  ranges  correspond 
to  a ray  path  oriented  at  an  elevation  angle  of  5°  and  intersecting  the  ionosphere  at  an  alti- 
tude of  approximately  100  and  300  km,  respectively.  For  this  configuration,  the  E-layer 
scattering  area  evaluates  to  a 95.  6 dBsm  and  79.  6 dBsm  for  a pulse  length  of  400  and  10  ns, 
respectively.  For  the  corresponding  pulse  lengths,  the  F-layer  scattering  area  is  99.  0 dBsm 
and  82.  9 dBsm.  It  is  seen  that,  except  for  the  upper  quartile  8. 125-MHz  results,  the  F- 
layer  cross  section  is  about  17  to  20  dB  greater  than  the  E -layer  at  8.  125  MHz  and  20  to 
23  dB  at  14.  875  MHz.  The  difference  between  the  F-layer  and  E-layer  cross  sections  are 
most  likely  due  to  the  fact  that  the  F-layer  scatterers  cover  a wider  spatial  extent  than  the 
E -layer  scatterers. 
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TABLE  3-13 

mr  BACK-SCATTER  RADAR  - ESTIMATED  AVERAGE  CROSS-SECTION 
OF  E-LA YeK  AND  F-LAYER  IRREGULARITIES  BASED  ON  POLAR 
FOX  II  DATA,  JANUARY  1972 


Frequency 

(MHz 

8.  125 


14.875 


Pulse 

Length 

0-13) 

400 


Average 
Cross  Soction 
(cJBsm) 


Statistical 

Parameter 

E- Layer 

F-Layer 

(dB) 
F- Layer /E 

Upper 

Quartile 

61.  2 

7 i • B 

t A 

* ><>  U 

Median 

48.  0 

65.4 

17.4 

Lower 

Quartile 

37.  0 

57.  0 

20.0 

Upper 

Quartile 

45.4 

55.7 

10.3 

Median 

32.  0 

40,  3 

17.  3 

Lower 

Quartile 

21.0 

40.9 

19.9 

Upper 

Quartile 

39.1 

61.  9 

22.8 

Median 

31.  7 

51.  8 

20.1 

Lower 

Quartile 

21.6 

44.7 

-23. 1 

Upper 

Quartile 

23.  1 

46.8 

22.7 

Median 

16.7 

35.7 

20.0 

Lower 

Quartile 

6.6 

28.6 

23.0 

Note-  Propagation  losses  not  removed  from  average  cross  section. 


prsT  AVAILABLE  COPY 


Figures  3-10  and  3-11  are  histograms  of  the  mean  cross  section  per  Hz  Doppler  band- 
width of  ionospheric  backscatter  echoes  deduced  from  Polar  Fox  II  measurements  made  at 
a frequency  of  8. 125  and  14.  875  MHz,  respectively,  in  July  1972  between  2000  and  0300 
hours  local  time.  The  data  which  were  supplied  by  the  Lincoln  Laboratory  for  incorpora- 
tion in  a report  by  Herman  and  Vargas-Vila  (1973)  are  based  on  the  assumption  that  the 
echoes  originated  from  a point  source.  It  should  be  noted  that  the  ionospheric  loss  term 
has  not  been  removed  from  the  calculations.  The  data  are  applicable  to  reflections  obtained 
in  the  azimuth  interval  from  -30°  to  7°  and  range  interval  from  624  to  2199  km.  In  addition, 
the  data  contain  various  modes  of  propagation  such  as  slant-F,  direct-F,  1-hop  F-layer 
plus  slant  and  direct  F,  direct-E  and  single-hop  E-layer  returns. 

The  theoretical  normal  distribution  curve  also  plotted  in  Figures  3-10  and  3-11  were 
derived  utilizing  the  statistical  parameters  of  the  experimental  data  as  specified  by  Herman 
and  Vargas-Vila  (1973).  At  8. 125  MHz,  the  mean  of  the  distribution,  n = 63. 3 dBsm/Hz  and 
the  standard  deviation,  a = 10.  0 dBsm/Hz  while,  at  14.  875  MHz,  \i  = 36.  7 dBsm/Hz  and 
cr  = 15.  0 dBsm/Hz. 

Evidence  that  the  8. 125-  and  14. 875-MHz  cross  section  data  are  approximately  log- 
normally  distributed  is  presented  in  Figures  3-12  and  3-13.  As  illustrated  in  Figure  3-4, 
a similar  distribution  is  obtained  for  the  amplitudes  of  the  401-MHz  E-layer  reflections  ob- 
served by  SRI  in  Scotland. 

According  to  Pollon  (1970),  no  satisfactory  physical  mechanism  for  the  occurrence 
of  log-normally  distributed  scattering  has  as  yet  been  found. 

The  statistical  parameters  of  the  radar  cross  section  listed  in  Table  3-14  are  ob- 
tained from  the  cumulative  distribution  functions  presented  in  Figures  3-12  and  3-13.  The 
difference  between  the  median  and  the  standard  deviation  deduced  by  this  analysis  and  by 
Herman  and  Vargas-Vila  (1973)  is  due  to  the  fact  that  the  exact  values  of  Herman's  and 
Vargas-Vila's  experimentally  derived  radar  cross  sections,  shown  in  Figures  3-10  and 
3-11,  were  not  available. 
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NO.  OF  OCCURRENCES 


NO.  OF  OCCURRENCES 


111516 


NORMAL  DISTRIBUTION 


70  80 


a/Lp  (dBsm/IIz) 


Figure  3-11.  Radar  Cross  Section  of  Polar  Fox  II  Ionospheric  Echoes 
Recorded  at  a Frequency  of  14.  875  MHz,  July  1972  (after 
Herman  and  Vargas-Vila,  1973) 
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TABLE  3-14 


itADAR  CROSS  SECTION  OF  FIELD-ALIGNED  IONIZATION  FOR  THE  AZIMUTH 
INTERVAL  FROM  -30°  TO  7°  AND  RANGE  INTERVAL  FROM 
624  TO  2199  km,  JULY  1972,  2000  TO  0300  HOURS  LOCAL  TIME 


Radar  Cross  Section  (dBsm/Hz) 

8. 125  MHz 

14.875  MHz 

Statistical 

Parameter 

GE  Analysis 

Herman  and 
Vargas-Vila 
Analysis 

GE  Analysis 

Herman  and 
Vargas-Vila 
Analysis 

Upper 

Decile 

74.5 

- 

56.  9 

- 

Upper 

Quartile 

69.2 

- 

47.3 

- 

Median 

62.3 

63.4 

34.6 

35.6 

Lower 

Quartile 

56.4 

- 

23.5 

- 

Lower 

Decile 

49.3 

- 

16.1 

- 

Standard 

Deviation 

9.8 

10.0 

16.8 

15.0 

Mean 

- 

63.  3 

- 

36.7 

Note:  Propagation  losses  not  removed  from  radar  cross  section. 
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CUMULATIVE  DISTRIBUTION  {%) 


a/Lp  (dBsm/Hz) 


Figure  3-12.  Cumulative  Distribution  Function  of  Radar  Cross  Section  of 
Polar  Fox  II  Ionospheric  Fchoes  Recorded  at  a Frequency  of 
8.  125  MHz,  July  1972 
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It  is  of  interest  to  note  the  independency  of  radar  cross  section  with  azimuthal 
orientation  and  range  was  found  to  exist  in  the  July  1972  ionospheric  echo  data  (Herman  and 
Vargas-Vila,  1973).  As  shown  in  Tables  3-10  and  3-11,  the  January  1972  data  also 
disclosed  the  azimuthal  independency. 

It  is  reiterated  that  the  Polar  Fox  II  estimates  of  the  E-  and  F-layer  scattering 
coefficients,  Tables  3-10  and  3-11,  and  radar  cross  sections,  Table  3-14, contain  the 
propagation  loss  terms,  i.  e. , the  ionospheric  propagation  losses  have  not  been  removed 
from  the  calculations.  Because  of  this,  it  is  not  possible  to  use  the  Polar  Fox  II  iono- 
spheric backscatter  echo  data  for  extrapolation  to  other  frequencies  in  the  HF  bu»jH. 


! 
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SECTION  IV 

DOPPLER  FREQUENCY  SPECTRUM  OF  IONOSPHERIC  BACKSCATTER  ECHOES 


The  Doppler  speetrum  of  radar  pulses  refleeted  from  the  aurora  differs  from  that  of 
the  original  pulse.  The  changes  that  ean  take  place  are:  (1)  the  center  frequency  can  be 
shifted;  and  (2)  the  spectral  width  can  be  increased. 

The  shift  in  the  center  frequency  corresponds  to  a drift  motion  of  the  auroral  ioniza- 
tion. Radar-auroral  data  indicate  an  east-west  drift  before  magnetic  midnight  and  a west- 
east  drift  after  magnetic  midnight  (Leadabrand,  et  al. , 1965).  Drift  velocities  of  the  order 
of  500  m/s  are  typical.  According  to  the  analysis  of  E-layer  radar-auroral  echoes  recorded 
in  Scotland  at  a frequency  of  401  MHz,  the  maximum  Doppler  shift  is  normally  ±2. 15  kHz 
(Larson  and  Hodges,  1967). 

The  spread  in  the  Doppler  spectrum  is  due  to  the  random,  turbulent  motion  of  the  ir- 
regularities of  electron  density  in  the  auroral  ionization.  Figure  4-1  which  is  the  Doppler 
frequency  speetrum  obtained  by  averaging  over  17  consecutive  pulses  (Leadabrand,  et  al. , 
1965)  gives  evidence  that  the  auroral  speetrum  is  significantly  broadened  compared  to  the 
transmitted  pulse  spectrum.  According  to  Larson  and  hedges  (1967),  the  maximum  Doppler 
frequency  spread  at  401  MHz  is  on  the  order  of  1.  90  kHz  at  -6  dB  below  the  peak,  4.  00  kHz 
at  -12  dB  and  4.  75  kHz  at  -18  dB. 

The  expected  (one  standard  deviation)  of  Doppler  frequency  shift  and  spread  at  HF 
are  presented  in  Figure  4-2.  The  estimates  are  based  on  the  fact  that  the  Doppier  fre- 
quency variations  are  directly  proportional  to  frequency  and  on  the  assumption  that  shifts 
and  spreads  are  normally  distributed.  Thus,  the  maximum  shift  and  spread  are  equivalent 
to  the  3-sigma  (standard  deviation)  value  of  the  distribution.  At  10  MHz,  the  Doppler  fre- 
quency shift  evaluates  to  ±18  Hz  and  ;he  Doppler  frequency  spread  becomes  16  Hz  at  the 
-6  dB  level,  33  Hz  at  -12  dB  and  40  Hz  at  -18  dB. 

The  estimates  of  E-layer  HF  radar-auroral  total  Doppler  frequency  deviationare 
plotted  in  Figure  4-3.  At  10  MHz,  it  is  possible  that  the  1-sigma  value  of  the  total  Doppler 
deviafion  eould  be  on  the  order  of  ±26  Hz  at  the  -6  dB  level,  ±35  Hz  at  -12  dB  and  ±38  Hz 
at  -18  dB. 
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Deviation  Based  on  SRI-Scotland  Measurements  at  401  MH 
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E-layer  radar -auroral  echoes,  having  radial  velocities  as  high  as  1450  m/s,  have 
been  observed  at  17  MHz  by  Brooks  (1906).  This  corresponds  to  a Doppler  frequency  shift 
(Ej)  of  104.  3 Hz  which  is  derived  from  the  expression 


where  V is  the  radial  velocity  and  X is  transmitted  wavelength. 

The  measurements  of  the  spectral  characteristics  of  F-layer  irregularities  made  under 
the  Polar  Fox  II  program  reveal  that  the  spectrum  is  somewhat  flat  over  the  ±15  Hz  Doppler 
bandwidth  at  frequencies  up  to  and  including  19.  25  MHz  (Campbell,  et  al. , 1972). 

The  maximum  radial  velocity  of  F-region  field-aligned  irregularities  detected  at 
17.  3 MHz  (Baggaley,  1970)  was  found  to  be  on  the  order  of  165  m/s  with  a mean  of  65  m/s. 

In  terms  of  the  Doppler  frequency  shift,  these  velocities  evaluate  to  19.  0 and  7.  5 Hz, 
respectively. 


SECTION  V 


FREQUENCY  OF  OCCURRENCE  OF  IONOSPHERIC  BACKSCATTER  ECHOES 
AND  CORRELATION  WITH  SOLAR-GEOPHYSICAL  CONDITIONS 

The  experimental  observations  of  field -aligned  echoes  at  19.  4 MHz  conducted  over  a 
three-year  period  (1961-19C3)  at  a site  loeated  in  the  vicinity  of  Boston  disclosed  that  the 
echoes  were  present  for  as  long  as  11  h/day  (Malik  and  Aarons,  1964).  Although  a seasonal 
pattern  was  elcarly  defined,  there  appeared  to  be  a tendency  for  the  reflections  to  oceur  on 
a greater  number  of  days  during  the  summer  than  during  the  winter.  The  correlation  of 
ceho  activity  with  sunspot  number  was  found  to  exist.  For  example,  in  1961  when  the  sun- 
spot number  was  53,  eehoes  appeared  on  73%  of  the  days  at  an  average  of  3.  0 h/day.  In 
1963  when  the  sunspot  number  was  28,  55%  of  the  echoes  appeared  at  an  average  of  1.  5 h/day. 
Maximum  ceho  activity  occurred  between  1800  and  2000  h local  time. 

A more  thorough  analysis  of  the  19.  4 MHz  field-aligned  echoes  recorded  over  a 
period  from  1961  through  1965  has  been  performed  by  the  Air  Force  Cambridge  Research 
Laboratories  (AFCRL)  (Basu  et  al. , 1973,  1974).  It  was  found  that,  under  quiet  magnetic 
conditions,  field-aligned  E-layer  eehoes  showed  a summer  evening  maximum  and  appeared 
to  be  associated  with  the  ground -backs  catter  echoes  from  sporadic  E (Basu,  et  al. , 1973), 

A weak  secondary  maximum  existed  in  the  winter  with  no  detectable  field -aligned  eehoes 
during  the  daytime.  During  disturbed  magnetic  conditions,  the  field-aligned  E-layer  eehoes 
increased  with  magnetic  activity  and  appeared  during  the  daytime  with  a corresponding 
decrease  in  sporadic  E.  The  data  did  not  reveal  a seasonal  dependence. 

Under  quiet  magnetic  conditions,  the  field-aligned  F-layer  echoes  were  found  to  be  a 
sunset  p U;»iorr.enon  and  correlated  with  the  sunspot  cycle  (Basu,  et  al. , 1974).  In  general, 
the  echoes  occurred  in  the  range  interval  from  1050  to  1500  km.  However,  during  the  day- 
time, long  range  eehoes,  at  times,  were  observed  at  ranges  of  3000  to  3300  km.  An  inter- 
esting disclosure  was  the  fact  that  the  occurrence  of  the  F-layer  field-aligned  echoes  in- 
creased directly  with  the  increase  in  magnetic  activity  until  the  level  Kpr  (Fredericksburg, 
Maryland,  three-hourly  K index)  attained  a value  of  four.  Beyond  this  level,  there  was  a 
decrease  in  the  occurrence  of  the  F-layer  echoes  with  no  echoes  being  observed  when 
Kpr  > 7.  Radar  auroral  measurements  made  at  Stanford  University  at  a frequency  of 
17. 3 MHz  reveal  that  maximum  auroral  activity  occurred  between  1900  and  0400  hours  loeal 
time  (Peterson  et  al. , 1955). 
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According  to  the  IIP  radar  echo  observations  conducted  at  12  and  18  MHz  at  lhillman, 
Washington,  there  was  a pronounced  peak  in  the  frequency  of  occurrence  of  the  echoes  from 
field-aligned  irregularities  near  the  time  of  local  sunset  (Hower,  et  al. , 19GG).  The  charac- 
teristics of  the  echoes  detected  at  sunset  were  different  from  those  occurring  at  night.  It 
was  found  that  the  local  sunset  echoes  appeared  at  ranges  of  about  2000  km  and  were  gen- 
erally spread  in  range  (diffuse).  The  occurrence  of  the  sunset  echoes  was  relatively  inde- 
pendent of  magnetic  activity.  The  nighttime  echoes,  on  the  other  hand,  appeared  at  ranges 
of  about  1000  km  and  tended  to  display  little  spread  in  range  (discrete).  The  occurrence  of 
nighttime  echoes  was  highly  dependent  on  magnetic  activity. 

Hower,  et  al. , (19GG)  noted  a decreasing  echo  activity  with  solar  cycle,  the  maximum 
percentage  occurrences  during  the  nighttime  being  approximately  80%  in  December  1958, 

37%  in  December  19G3  and  35%  in  March  19G4.  The  (Zurich)  relative  sunspot  number  during 
the  corresponding  period  was  185.  2,  11.  8 and  14.  5,  respectively. 

Hower  and  Makhijani  (19G9)  have  concluded  that  HF  F-layer  baekscatter  echoes  arise 
from  the  same  general  Irregularity  regions  in  which  spread- F is  detected. 

Sprenger  and  Glode  (19G4)  have  reported  that  the  diurnal  variation  of  the  frequency 
of  occurrence  of  radar -auroral  echoes  recorded  from  late  1958  to  mid  1962  at  33  MHz  at 
Kuhlungsborn  (geographic  coordinates:  54°N,  12’E;  geomagnetic  latitude:  54°N)  showed  a 
double-maximum  at  approximately  0100  and  1700  hours  local  time.  The  probability  of  oc- 
currence of  E -layer  auroral  echoes  was  found  to  be  a function  of  magnetic  activity.  No 
auroral  echoes  were  observed  at  a magnetic  index,  Kp  < 5 and  100%  occurrence  at  Kp  = 9. 
With  regard  to  the  correlation  of  auroral  echo  activity  with  sunspot  cycle,  Sprenger  and 
Glode  (1964)  have  concluded  that  the  maximum  of  auroral  activity  occurred  in  1960  about 
two  years  after  the  sunspot  maximum  and  that  the  activity  decreased  to  zero  within  a period 
of  only  one  and  a half  years. 

Brooks'  (19G5,  19G6)  investigation  of  radar-auroral  echoes  at  17  MHz  conducted  in  the 
1959-1960  period  indicated  that,  for  the  discrete  echoes,  the  maximum  number  appeared 
between  0000-0300  hours  local  time  and,  for  the  diffuse  echoes,  between  1900-2100  hours 
local  time. 
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A maximum  of  occurrence  of  baekscattcr  echoes  at  13.  8GG  MHz  from  the  F-rrgion 
on  75%  of  the  nights  for  March  1958  was  reported  by  Weaver  (19G5).  A minimum  number 
of  echoes  were  observed  in  June,  the  data  being  collected  during  1957  and  1958  at  Ithaea, 
New  York,  with  the  antenna  beam  oriented  in  a northward  direction. 

Baekscattcr  observations  of  F-region  field-aligned  irregularities  made  at  17  MHz 
near  Sheffield,  England  (geographic  coordinates:  53.  43°N,  1. 58°W;  geomagnetic  latitude: 

5G.  4°N)  between  mid-October,  19G4,  and  mid -January,  19GG,  (Baggaley,  1970)  which  was  a 
period  of  low  sunspot  activity,  showed  that,  on  the  average,  the  echoes  were  present  on 
23.3%  of  the  days  while  only  1.  5%  of  the  total  observing  time  was  occupied  by  the  echoes. 

A summer  maximum  and  a winter  minimum  were  found  to  exist.  This  seasonal  variation 
correlates  with  Malik's  and  Aarons'  (19G4)  results. 

With  regard  to  the  temporal  variation,  Baggaley  (1970)  noted  that  there  was  a correla- 
tion between  the  onset  of  F-layer  echoes  and  the  solar  zenith  angle  at  the  reflection  point  in 
the  F-layer  (assuming  a height  of  300  km).  That  is,  a maximum  number  of  echoes  occurred 
at  a solar  zenith  angle  of  90-95°,  and  a minimum  at  angles  less  than  85°  and  greater  than 
125°.  In  addition,  no  correlation  was  found  to  exist  between  F-layer  echo  occurrence  and 
magnetic  activity.  Of  the  total  time  for  which  F-layer  echo  activity  was  present,  only  11% 
was  associated  with  magnetic  index  Kp  > 4. 

Field-aligned  echoes  from  the  F-region  observed  at  Brisbane,  Australia,  at  16  MHz 
occurred  preferentially  during  geomagnetic  disturbances  and  correlated  strongly  with 
spread-F  and  radio  star  scintillations  occurring  in  the  same  region  of  the  ionosphere 
(Swenson,  1972). 

An  analysis  of  29-MHz  baekscattcr  measurements  carried  out  in  Northern  (West) 
Germany  and  Seandanavia  (geomagnetic  latitudes  of  55-77°N)  by  Czechowsky,  et  al.  (1974) 
confirmed  the  close  correlation  between  radio-auroral  occurrence  and  geomagnetic 
activity.  Most  of  the  radio  auroras  appeared  in  the  afternoon  and  evening  hours. 

The  Polar  Fox  II  HF -ionospheric  baekscattcr  measurements  revealed  that,  in  the 
summertime,  the  echoes  appeared  most  often  at  night  and  at  the  shorter  ranges.  During 
the  wintertime,  the  long  range  echoes  appeared  most  often  during  the  daytime  while  the 
short  range  echoes  were  equally  present  during  the  day  and  night.  The  maximum  number  of 
echoes  appeared  in  the  northerly  direction.  The  echoes  would  appear  at  a given  range  at 
one  frequency  and  not  at  another  (Edwards,  et  al.  , 1973). 

Table  5-1  summarizes  the  experimental  results  of  the  characteristics  of  HF  iono- 
spheric baekseatter  echoes  and  their  correlation  with  solar-geophysical  conditions. 
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Loral  Time 


SECTION  VI 


CON  C I .US  IONS 


The  spatial  distribution  of  II F radar-ionospheric  backseatter  echoes  can  be  readily 
determined  utilizing  the  concept  that  the  reflections  originate  in  regions  of  space  where  the 
direction  of  propagation  is  near  perpendicular  to  the  direction  of  the  lines  of  force  of  the 
earth's  magnetic  field,  and  where  high  auroral  activity  and  ionospheric  irregularities  are 
encountered. 

It  is  estimated  that,  for  a midlatitude  location,  II F E-layer  auroral  echoes  should 
be  nonexistent  during  the  daytime  while,  during  the  evening  for  quiet  magnetic  conditions 
(Q  - 1),  E -layer  echoes  could  appear  between  the  azimuths  of  -39c  and  +19°  and  between 
the  ranges  of  5fi0  and  780  nmi.  Under  severe  magnetic  conditions  (Q  = 7),  the  azimuthal 
and  range  extents  could  increase  to  -02  to  +20c  and  150  to  780  nmi,  respectively. 

E-layer  echoes  can  be  expected  to  be  observed,  during  both  the  daytime  and  nighttime, 
out  to  ranges  on  the  order  of  1250  nmi,  while  the  azimuthal  extent  of  the  F-layer  reflections 
could  cover  the  entire  region  between  -90J  to  +903. 

In  the  frequency  range  between  5-30  MHz,  the  median  E-layer  radar-auroral  SNR 
deduced  from  UIIF  radar-auroral  data  acquired  by  SRI  at  Scotland,  could  vary  between 
88-94  dB  and  72-78  dB  for  a pulse  length  of  400  and  10  ps,  respectively,  with  the  maximum 
occurring  at  10  MIIz.  The  quartile  values  of  the  E-layer  SNR  arc  ±10  dB  with  respect  to  the 
median  level  while  the  decile  values  arc  ±19  dB. 

An  analysis  of  the  IIF  ionospheric  backseatter  data  obtained  in  Maine  under  the 
Polar  Fox  II  program  reveals  that,  except  for  the  upper  quartile  8.  125-MIIz  results,  the 
F-layer  cross  section  could  be  17-20  dB  greater  than  the  E-layer  cross  section  at  8. 125 
MIIz  and  20-23  dB  at  14.  875  MHz. 

It  is  found  that  the  amplitude  of  E-layer  auroral  echoes  and  the  cross  section  of  E- 
and  F-layer  reflections  arc  log-normally  distributed. 


Kxtrapolating  from  the  SRI-Scotland  UHF  radar-auroral  data,  it  found  that,  for 
5 MHz  radar  transmissions  reflected  from  K-layer  auroral  ionization,  the  1-sigma  value 
of  the  total  Doppler  frequency  deviation,  i.  e. , Doppler  frequency  shift  and  spread,  could 
be  on  the  order  of  ±13  Hz  at  the  -6  dB  normalized  signal  level,  ±17  Hz  at  -12  dB  and  ±19  Hz 
at  -IS  dB,  At  30  MHz,  the  1-sigma  Doppler  deviation  should  increase  to  approximately 
±78  Hz,  ±105  Hz  and  ±114  Hz,  respectively. 

Based  on  the  Polar  Fox  II  HF  ionospheric  measurements,  the  Doppler  frequency 
spectrum  of  field-aligned  F-layer  backscatter  appears  to  be  flat  over  the  ±15  Hz  bandwidth 
up  to  a frequency  of  19.  25  MHz. 

The  frequency  of  occurrence  of  HF  field-aligned  backscatter  echoes  maximizes,  for 
the  most  part,  near  local  sunset.  However,  experimental  observations  have  shown,  at 
times,  peak  activity  near  the  midnight  hours. 

There  is  a strong  correlation  of  backscatter  echoes  with  the  solar  cycle.  That  is, 
the  percentage  of  days  displaying  backscatter  reflections  decrease  with  decreasing  sunspot 
number. 

The  characteristics  of  sunset  echoes  are  found  to  be  different  from  those  occurring 
during  the  night.  The  sunset  echoes  are  generally  of  the  diffuse  type  and  are  Independent 
of  magnetic  activity  which  the  nighttime  echoes  are  discrete  and  correlate  with  magnetic 
activity. 
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APPENDIX  A 


RADAR-AURORAL  CLUTTER  (RAC)  COMPUTER  PROGRAM 

The  radar-auroral  clutter  (RAC)  computer  program  provides  a means  for  determining: 
(1)  the  magnetic  field  geometry  necessary  for  radar-auroral  backscatter  reflections  to  take 
place,  i.  e.  , angular  orientation  (aspect  angle)  of  the  direction  of  propagation  with  respect 
to  the  direction  of  the  lines  of  force  of  the  earth's  magnetic  field,  and  (2)  the  locations  on  the 
earth's  surface  for  the  reception  of  magnetic  field-forward  scatter  from  a given  transmitter 
location. 

The  detailed  mathematical  developments  of  the  backscatter  and  forward  scatter  mag- 
netic field  geometries  have  been  published;  thus  they  will  not  be  discussed  in  this  report 
(Millman,  1958a,  1959,  1959).  In  the  case  of  radar-auroral  backscatter,  the  magnetic 
field-aspect  angle,  0,  can  be  expressed  by  the  function 

0 - f(<>,  X,  A,  E,  h,  I,  D,  nt,  N’e> 

where  <>  and  X are  the  geographic  latitude  and  longitude  of  the  radar  site,  respectively,  A 
and  E are  the  azimuth  and  elevation  angles  of  the  antenna  beam,  respectively,  h is  the 
altitude  above  the  earth's  surface  of  the  reflection  point,  I and  D are  the  magnetic  Indication 
and  declination  angles,  respectively,  n is  the  index  of  refraction  in  the  troposphere  and  N 

l t 

is  the  electron  density  in  the  ionosphere.  A more  complete  discussion  of  this  function  is 
given  in  Appendix  B. 

The  RAC  computer  program  has  the  capability  of  determining  the  orientation  of  the 
lines  of  force  of  the  earth's  magnetic  field  (angle  0)  for  any  geographic  location,  antenna 
azimuth  angle,  elevation  angle  and  radar  range  or  altitude  above  the  earth's  surface. 

An  option  is  provided  in  the  program  to  employ  one  of  the  following  three  models  of 
the  earth's  magnetic  field: 

1.  Graphical  model 

2.  Dipole  model 

3.  Spherical  harmonic  model 


A-l 


In  the  graphical  model,  the  parameters  I and  I)  have  been  sealed  in  2.  52 * *  5 and  5°  latitude 
and  longitude  increments  over  the  whole  earth’s  surface  from  the  epoch  1955  isogonic  and 
isoclinic  maps  issued  by  the  U.  S.  Navy  Hydrographic  Office  and  the  Canadian  Department  of 
Mines  and  Technical  Surveys.  The  ground  observed  magnetic  data  obtained  from  these  maps 
are  stored  in  matrix  form  in  the  11(135  digital  computer.  Linear  interpolation  is  used  for 
determining  the  values  of  the  magnetic  field  elements  at  other  geographic  locations.  For 
propagation  through  the  ionosphere,  the  angles  I and  D at  the  subionospheric  point  are  as- 
sumed to  be  invariant  at  all  altitudes.  According  to  the  1955  magnetic  charts,  the  geo- 
graphic positions  of  the  north  and  south  dip-poles,  which  are  dependent  upon  local  surface 
anomalies  and  irregularities,  are  approximately  73.  5’N’,  101.  1°\V  and  67.  6®S,  144.  0°K, 
respectively. 

In  the  dipole  model,  it  is  assumed  that  the  earth's  magnetic  field  can  be  approximated 
by  a magnetic  dipole  located  at  the  earth's  center  having  geomagnetic  poles  at  78.  6°N, 

70.  1°W,  and  78.  6°S,  109.  9°E  (Vestine,  et  al. , 1947).  For  the  magnetic  dipole  approxima- 
tion, the  inclination  and  declination  angles  are  invariant  with  altitude. 

The  spherical  harmonic  model  assumes  that  the  earth's  main  field,  i.  e. , the  mag- 
netic field  that  excludes  such  phenomena  as  magnetic  disturbances  and  diurnal  variations, 
can  be  represented  by  a regular  or  dipole  field  and  an  irregular  field.  A series  of  spherical 
harmonies  expresses  the  magnetic  potential  function  for  the  main  field  over  the  earth's 
surface,  the  spherical  harmonic  coefficients  being  input  parameters  in  the  RAC  program. 

For  the  spherical  harmonic  model,  the  angles  I and  D are  altitude  dependent.  The  mathe- 
matical formulation  of  the  dipole  and  spherical  harmonic  models  are  contained  in  Appendix  B. 

The  calculation  of  the  angle  0 is  performed  for  either 

1.  An  unrefraeted  ray  which  is  applicable  for  line-of-sight  radars  operating  at  a 
frequency  of  approximately  400  MHz  and  above  (Millman,  1972,  1975)  and  with 
the  antenna  beam  elevated  at  angles  greater  than  about  5°  (Millman,  1975)  or 

2.  A refracted  ray  which  takes  into  account  the  deflection  introduced  by  the  tropo- 

sphere and/or  ionosphere. 

The  refractive  bending  in  the  troposphere  is  small  compared  to  the  ionospherie  effect  at  HF 

(Millman,  1975)  and  thus,  as  a first  approximation,  is  neglected  in  the  magnetic  field  com- 
putations appropriate  to  HF  radars. 


The  basic  assumption  used  in  the  refraction  computation  is  that  the  troposphere  and 
ionosphere  can  be  considered  to  be  stratified  into  spherical  layers  of  constant  refractive 
index.  The  computational  procedure  has  already  appeared  in  the  literature  (Millman,  1958b) 
and  therefore  will  not  be  discussed  in  this  report.  The  accuracy  of  this  method,  which  is 
dependent  upon  the  number  of  stratified  layers,  can  be  greatly  improved  by  merely  in- 
creasing the  numbers  of  stratified  layer,  i.  e. , decreasing  the  thickness  of  each  individual 
layer  element  (Millman,  1951). 

The  computer  program  incorporates  the  option  for  single  or  multiple  ionospheric 
"hops",  i.  e. , ionospheric  reflection  points,  with  a different  ionospherie  eleetron  density 
profile  for  each  individual  "hop".  The  magnetic  field  aspect  angle  is  computed  along  both 
the  upward  ray  and  the  downward  ray  at  every  "hop". 

An  option  is  available  in  the  program  for  determining,  in  the  unrefraeted  ray  ease, 
the  azimuth  angle  corresponding  to  a specified  preset  aspect  angle  ft  for  a set  of  values  of 
elevation  angles  and  radar  ranges  or  ionospheric  altitudes.  For  the  refracted  ease,  this 
option  provides  a means  for  obtaining,  along  any  ray  path,  the  altitude  at  whieh  a specified 
aspect  angle  is  attained. 

The  probability  of  observing  radar  reflections  from  the  aurora  is,  in  addition  to  the 
geometric  and  height  restriction,  dependent  upon  the  frequency  of  occurrence  of  auroral 
activity  in  the  region  of  interest.  Auroral  reflections  can  be  expected  to  exist  in  regions 
where  both  the  conditions  of  near-perpendieularity  at  ionospherie  heights  and  high  auroral 
activity  are  satisfied.  The  fulfillment  of  only  one  requirement  is  not  sufficient  to  warrant 
a radar-auroral  return.  The  auroral  activity  model  of  Vestine  (1944),  shown  in  Figure  A-l, 
is  presently  an  integral  part  of  the  RAC  program.  The  contour  plots,  whieh  are  stored  in 
matrix  form  in  the  program,  give  the  estimated  frequency  of  days  with  occurrence  of  visual 
aurora,  on  clear,  dark  nights  in  the  northern  hemisphere.  The  probability  of  auroral  oc- 
currence is  computed  along  the  ray  paths  for  both  the  unrefraeted  and  refracted  rays.  It  is 
realized  that  Vcstinc's  model  is  quite  restricted  since  it  applies  only  to  the  midnight  hours 
and  docs  not  take  into  account  the  effect  of  magnetic  activity.  A more  appropriate  model  to 
use  for  this  application  is  the  Feldstein’s  and  Starkov’s  (1967)  auroral  belt  model,  presented 
in  Figure  A-2,  mainly  because  the  spatial  orientation  of  the  auroral  activity  region  is  both 
time-  and  geomagnetic  activity  dependent. 

Figure  A-3  is  a logic  diagram  of  the  RAC  computer  program.  The  dashed  line  around 
the  Feldstein-Starkov  auroral  belt  model  signifies  that  this  function  is  not  presently  incor- 
porated in  the  RAC  program. 
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Figure  A -3.  Radar -Auroral  Clutter  (RAC)  Computer  Program  Logic  Diagram 


APPENDIX  B 

MAGNETIC  FIE LD-BACKSCATTER  GEOMETRY 


In  the  study  of  radar-auroral  reflections,  it  is  desirable  to  determine  the  orientation 
of  the  lines  of  force  of  the  earth’s  magnetic  field  for  any  geographic  location,  antenna 
azimuth  bearing,  elevation  angle,  and  reflection  height  above  the  earth's  surface. 


It  can  be  shown  that  the  propagation  angle  0 (i.  e.  , the  angle  between  th»  direction  of 
the  magnetic  lines  of  force  and  the  direction  of  electromagnetic  propagation),  depicted  in 
Figure  B-l,  can  be  expressed  by  the  relationship  (Millman,  1959; 

0 = cos  * (-cos  e sin  I - sin  e cos  I cos  (y  - D)  ] (B-l) 


where  I and  D are  the  magnetic  inclination  and  declination  angles,  respectively.  These 
parameters  specify  the  direction  of  the  total  magnetic  intensity  veetor  in  space. 

The  angle  y whieh  is  the  geographic  azimuth  bearing  of  the  observer's  location  as 
measured  at  the  subionospheric  point,  i.  e.  , the  location  on  the  earth‘s  surface  directly  be- 
neath the  magnetic  field  intersection  point,  is  defined  by  (Millman,  195^a) 

. f sin  (X  - Xp)  cos  4>n 

y = tan"  -7— r t — — r r r — r (B— 2) 

' sin  4>r  cos  4>p  - cos  <J>R  sin  <*>p  eos  (XR  - Xp)  v ' 

where  4>  and  X are  the  geographic  latitude  and  east  longitude,  respectively.  The  subscripts, 
R and  P,  refer  to  the  transmission  (or  radar)  site  and  the  reflection  point,  respectively. 


The  angle  e,  whieh  is  the  angle  between  the  ray  path  and  the  zenith  at  the  point  of 
reflection,  is  given  by 


e 


o 

+ h 


eos  E 


(B— 3) 


where  rQ  is  the  radius  of  the  earth,  E is  the  elevation  angle  of  the  antenna  beam  and  h is  the 
altitude  above  the  earth's  surface. 
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Figure  B-l.  Radar-Auroral  Geometric  Configuration 


For  the  dipole  model  (i.  c. , earth's  magnetic  field  approximated  by  a magnetic  dipole 
having  the  geomagnetic  north  pole  at  78.  6°N,  70. 1°W)  the  inclination,  or  magnetic  dip,  which 
specifies  the  direction  of  the  total  magnetic  intensity  vector  with  respect  to  the  horizon  and 
is  measured  positive  in  the  downward  direction,  is  given  by  (Chapman,  1951). 

1 - tan-1  [2  tan  tip]  (B-4) 

where  tl  is  the  geomagnetic  latitude. 

A simple  transformation  for  expressing  t/  p in  terms  of  geographic  coordinates  is  given 
by  (Millman,  1959) 

tip  ~ sin-1  (cos  (XM  - Xp)  eos  4>^.  cos  <t>p  + sin  sin  4>p]  (B-5) 

where  the  subscript,  M,  refers  to  the  geomagnetic  pole. 

The  magnetic  declination  angle,  which  specifies  the  direction  of  the  geomagnetic  pole 
and  is  measured  positive  in  a clockwise  direction  from  geographic  north,  can  be  defined  by 
the  function  (Millman,  1958a) 


sin  (XM  - Xp)  cos  4>m 


cos  4>p  sin  <$>M  - cos  (XM  - Xp)  sin  4>p  cos  <J>j 


It  is  of  interest  to  note  that,  for  the  dipole  model  approximation,  the  inclination  and 
declination  are  invariant  with  altitude. 

Due  to  local  magnetic  anomalies  on  the  earth's  surface,  there  arc  large  discrepancies 
and  widespread  departures  from  a simple  dipole  field.  To  take  this  into  account,  the  earth's 
main  field,  i.  e. , the  magnetic  field  which  excludes  such  phenomena  as  magnetic  disturbances 
and  diurnal  variations,  can  be  represented  by  a regular  or  dipole  field  and  an  irregular  field. 
A magnetic  potential  function  for  the  main  field  over  the  earth's  surface  can  be  expressed  in 
terms  of  a scries  of  spherical  harmonics  (Chapman  and  Bartels,  1940)  such  as 

n = « m-n  . . n + 1 

V.  £ ^2  (gnm  COS  (mX) + hnm  sin  (mX)]  P^(cos  4>')  (B-7) 

n = 0 m = 0 

where  r is  the  distance  from  the  center  of  the  earth,  4>'  is  the  geographic  colatitude, 

Pm  (cos  <>')  is  the  associated  Legendre  function  of  degree  n and  order  m and  g and  h 
are  the  coefficients  of  the  spherical  harmonic  expansion. 
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It  should  be  noted  that  the  first-degree  harmonic  terms  in  Equation  (11-7),  i.e.  , terms 
with  n 1 and  :r>  = 0,  reduce  to  that  of  a dipole  potential. 

In  practice,  the  magnetic  potential  is  not  a measurable  quantity.  The  quantities  which 
are  measured  arc  namely  X,  the  northward  horizontal  component,  Y,  the  eastward  horizon- 
tal component,  and  Z,  the  downward  vertical  component  of  the  total  magnetic  intensity, 
which  are  defined  in  terms  of  the  magnetic  potential  by 

1 3 V 

x " F W (B-8) 


r sin  <$>'  3X 


y _ av 

L dr 


(B-10) 


By  definition,  the  magnetic  inclination  and  declination  angles  are  related  to  components  of 
the  total  magnetic  intensity  by 


I = tan" 


2 2 2 
(X  + Y ) 


(B-Il) 


D = tan 


"M 


(B-12) 


It  should  be  evident  from  Equations  (B-7)  through  (B-12)  that,  for  the  spherical  harmonic- 
magnetic  field  model,  the  parameters  I and  D are  altitude  dependent. 

For  a given  set  of  spherical  harmonic  coefficients,  the  magnetic  potential  function, 
Equation  (B-7),  becomes  specifiable.  This  function,  in  turn,  is  then  employed  to  evaluate 
the  spatial  distribution  of  I and  D. 

The  spherical  harmonic  coefficient  models  that  are  incorporated  In  the  RAC  computer 
program  are  listed  in  Table  B-l, 
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TABLE  B-l 

SPHERICAL  HARMONIC  MODELS  INCORPORATED 
IN  THE  RAC  COMPUTER  PROGRAM 


Coefficients 

Model 

Type 

Normalization 

Epoch 

Degree 

n 

Order 

m 

Total 

No. 

CalT  and  Sweeney  (1970) 

Schmidt 

1968.0 

10 

10 

120 

Eouguerc  (1965) 

Gaussian  and 
Schmidt 

1960.0 

8 

7 

77 

Hendricks  and  Cain  (1966) 

Schmidt 

1960.0 

9 

9 

96 

Hurwitz  et  al.  (1966) 

Schmidt 

1965.0 

12 

12 

155 

International 
Geomagnetic  Reference 
Field  (1969) 

Schmidt 

1965.0 

8 

8 

79 

Jensen  and  Cain  (1962) 

Gaussian 

1960.  0 

6 

6 

48 

In  the  computation  of  the  magnetic  field-aspect  angle,  defined  by  Equation  (B-l),  48 
spherical  harmonic  coefficients  (n  = G,  m = G)  are  only  used  although  the  International  geo- 
magnetic reference  field  (IGRF)  derived  by  the  International  Association  of  Geomagnetism 
and  Aeronomy  (IAGA)  and  the  models  proposed  by  Cain  and  Sweeney  (1970),  Fouguere  (19G5), 
Hendricks  and  Cain  (19GG)  and  Hurwitz,  et  al.  (19GG)  consist  of  a greater  number  of  coef- 
ficients. The  error  introduced  by  neglecting  the  higher  degree  and  order  coefficients  is  of 
minor  importance  in  geometric  studies  of  field-aligned  ionization. 

According  to  Cain  and  Sweeney  (1970)  the  IGRF  model  is  inadequate  when  a high  de- 
gree of  accuracy  is  required  especially  at  altitudes  below  1000  km.  For  this  condition, 

Cain  and  Sweeney  (1970)  recommend  that  models  such  as  the  one  they  propose  be  used. 
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